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ABSTRACT 

We present results of a 100 ks XMM-Newton observation of the Seyfert 1.5 AGN NGC 3227. Our 
best-fit broadband model to the EPIC pn spectrum consists of a moderately flat (photon index of 1.57) 
hard X-ray power-law absorbed by cold gas with a column density of 3 x 10^^ cm^, plus a strong soft 
excess, modeled as a steep power law with a photon index of 3.35, absorbed by cold gas with a column 
density of 9 x 10^° cm~^. The soft excess increases in normalization by ~20% in ^^20 ks, independently 
of the hard X-ray emission component, and the UV continuum, tracked via the OM, also shows a 
strong increasing trend over the observation, consistent with reprocessing of soft X-ray emission. Warm 
absorber signatures are evident in both the EPIC and RGS spectra; we model two absorbing layers, 
with ionization parameters log^ ~ 1.2 and 2.9 erg cm s~^, and with similar column densities (~ 1-2 
xlO^^ cm~^). The outflow velocities relative to systemic of the high- and low-ionization absorbers are 
estimated to be -(2060li7°) s"^ and -(420+i^°) km s~\ respectively. The Fe Ka fine FWHM width 
is 7000 ± 1500 km s^^; its inferred distance from the black hole is consistent with the BLR and with the 
inner radius of the dust reverberation-mapped by Suganuma et al. An emission feature near 6.0 keV is 
modeled equally well as a narrow redshifted Fe K line, possibly associated with a disk "hot-spot," or as 
the red wing to a relativistically broadened Fe line profile. Swift-BAT and archival RXTE data suggest 
at most a weak Compton refiection hump {R < 0.5), and a high-energy cutoff near 100 keV. From RXTE 
monitoring, we find tentative evidence for a significant fraction of the Fe line flux to track variations in 
the continuum on time scales < 700 days. 

Subject headings: galaxies: active — galaxies: Seyfert — X-rays: galaxies — galaxies: individual (NGC 
3227) 



1. INTRODUCTION 

A 1993 Advanced Satellite for Cosmology and Astro- 
physics (ASCA) observation of the nucleus of the Seyfert 
1.5 NGC 3227 first revealed evidence for ionized absorb- 
ing gas with iVn.wA ~ 1-4x10^^ cm~^ (Netzer et al. 1994, 
Ptak et al. 1994). Komossa & Fink (1997b), using 1993 
Rontgen Satellite (ROSAT)-VSVG data, and George et al. 
(1998b), using ASCA data taken in 1993 and 1995, fur- 
ther modeled the warm absorber, finding it to be rela- 
tively lowly-ionized relative to most other Seyferts' warm 
absorbers. Given the low-resolution spectra, the solution 
for further progress was to obtain a gratings spectrum. 
Komossa (2002) reported on a C/iandra-LETGS observa- 



tion of NGC 3227 in October 2000, noting edge-hke fea- 
tures near 0.7 keV. XMM-Newton first observed NGC 3227 
in November 2000, when the source was undergoing a 3- 
month long period of very high levels of obscuration due 
to a compact cloud with A'h near 3 x 10^"^ cm~^, cover- 
ing 90% of the central source, and likely associated with 
the BLR, traversing the line of sight (Lamer et al. 2003). 
The soft X-ray data from this observation indicated the 
warm absorber to have ionization parameter log^3 =1-7- 
2.0 erg cm s~^ and a column density A'h.wa — 2-9x10^^ 
cm~2 (Gondoin et al. 2003). 

Optical spectra of the nucleus of NGC 3227 have indi- 
cated strong reddening due to the presence of dust along 



The ionization parameter ^ = Lion'^e ^, where Lion is the isotropic 1-1000 Ryd ionizing continuum luminosity, ne is the electron number 
density, and r is the distance from the central continuum source to the absorbing gas. 
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the line of sight. The narrow hne Hq;/H/3 ratio, for in- 
stance, indicates reddening, with measured ratios in NGC 
3227 implying visual extinctions Ay — (Hq!/H/3)/3.1 near 
1.2-1.7 (Cohen 1983; Gonzalez Delgado & Perez 1997), 
though values of 4.5-4.9 have been reported (Mundell 
et al. 1995a, Rubin & Ford 19680 The broad Ha/H/3 
line ratio shows a similar level of reddening. Ay 1.4 
(Cohen 1983, Winge et al. 1995). The steep drop in NGC 
3227's continuum flux from the near-UV to the far-UV also 
supports reddening: the spectral index from International 
Ultraviolet Explorer (lUE) measurements is aiuE = —2.9 
for NGC 3227, whereas aiuE = —1.4 is more typical for 
Seyfert Is (Komossa & Fink 1997b; Courvoisier & Paltani 
1992; Kinney et al. 1991). Assuming a standard Galactic 
gas/dust ratio of A^h/^v = 2x10^^ cm~^ magn~^ (e.g., 
ShuU & van Steenburg 1985), the amount of NLR and 
BLR reddening implies that the dust is associated with a 
column of gas with iVn 2-3 xlO^^ cm~^. 

Komossa & Fink (1997b) and George et al. (1998b) 
also reported local absorption due to cold material with 
a column density A^h = 3-6 xlO^° cm~^, in addition to 
the Galactic column (2.1x10^° cm~^. Dickey & Lockman 
1990, Murphy et al. 1996). The inferred local cold column 
density agrees with that derived from 21cm VLA observa- 
tions (Mundell et al. 1995b); however, it is not sufficient to 
produce the observed optical reddening, assuming a stan- 
dard Galactic gas/dust ratio. The agreement between the 
inferred column density associated with the dust and the 
column density of the ionized gas, however, prompted Ko- 
mossa & Fink (1997b) to suggest that the dust was em- 
bedded in the warm absorber and not associated with the 
cold gas. Brandt, Fabian, & Pounds (1996) first suggested 
this "dusty warm absorber" concept for the quasar IRAS 
13349-1-2438. 

However, Kraemer et al. (2000) suggested that it was 
unlikely dust could survive within the high-ionization gas. 
They proposed a configuration consisting of two warm ab- 
sorber zones, in addition to the cold column. One warm 
absorber is dust-free and highly-ionized, responsible for 
the oxygen edges; the other is a very lowly-ionized ab- 
sorber with dust mixed in, referred to as a "dusty luke- 
warm absorber." Those authors modeled the latter ab- 
sorber with a density of 20 cm~'^ and a location encom- 
passing or lying just outside the NLR, which has a radius 
of ~100 pc (based on [O III] imaging; Schmitt & Kinney 
1996). Using Huhhle Space Telescope (HST)-STIS, Cren- 
shaw et al. (2001) first observed the features of the dusty 
lukewarm absorber, evident in the form of optical/UV ab- 
sorption lines due to intermediate species of C, N and Si; 
Crenshaw & Kraemer (2001) suggested that the absorber 
may possibly be a relatively highly-ionized component of 
the NLR seen in absorption. Kraemer et al. (2000) sug- 
gested that the dust may evaporate off the putative molec- 
ular torus and be swept up in an outflowing wind (e.g., 
Reynolds 1997). Based on V-band-to-K-band reverbera- 
tion mapping (the latter band dominated by thermal dust 
emission) conducted by the MAGNUM collaboration, Sug- 
anuma et al. (2006) concluded that the innermost extent 
the dust was at a radius of '^5-20 light-days, suggesting 
an alternate or additional site from which an outflowing 



wind could sweep up dust. 

In 2006, we obtained a new, 100 ks XMM-Newton ob- 
servation of NGC 3227 in an unobscured state, to better 
characterize the broadband X-ray continuum and X-ray 
warm absorber features, search for signatures of dust in the 
X-ray spectrum, relate the X-ray and dusty optical/UV 
warm absorbers, and study Fe K bandpass emission fea- 
tures. §2 describes the observations and data reduction. 
§3 and §4 describe EPIC spectral fits to the time-averaged 
spectrum, focusing on the Fe K bandpass and then extend- 
ing to the 0.2-10 keV bandpass, respectively. §5 describes 
the fits to the RGS spectra. §6 describes spectral analysis 
of archival RXTE data to constrain the hard X-ray contin- 
uum. §7 describes time- and flux-resolved spectral fits to 
the EPIC data, including Fvar spectra, and time-resolved 
spectral fits to the RXTE data to quantify the temporal 
behavior of the Fe line. X-ray/UV correlations are pre- 
sented in §8. The results are discussed in §9, and a brief 
summary is given in §10. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. XMM-Newton data reduction 

NGC 3227 was observed by XMM-Newton during Rev- 
olution 1279, on December 3-4, 2006, for a duration of 
99 ks. XMM-Newton^ s European Photon Imaging Cam- 
era (EPIC) consists of one pn CCD back-illmninated ar- 
ray sensitive to 0.15-15 keV photons (Striider et al. 2001), 
and two MOS CCD front-illuminated arrays sensitive to 
0.15-12 keV photons (MOS 1 and MOS 2, Turner et al. 
2001). Data from the pn were taken in Large Window 
mode; data from both MOSes were taken in PrimePar- 
tialW2/small window mode. The medium filter was used 
for all three EPIC instruments. Spectra were extracted 
using XMM-Newton-S AS version 7.1.0 and using standard 
extraction procedures. For all three cameras, source data 
were extracted from a circular region of radius 40"; back- 
grounds were extracted from circles of identical size, cen- 
tered ~3' away from the core. Hot, fiickering, or bad pixels 
were excluded. Data were selected using event patterns 0- 
4 for the pn and 0-12 for the MOSes. Using the SAS task 
EPATPLOT, we verified that pile-up was negligible for the 
pn. For the MOS, pile-up was minimal (<3%) up to 8 
keV, and a bit higher (3-10%) above 8 keV. The effect of 
such pile-up is to artifically fiatten the spectrum. How- 
ever, we also extracted counts from the MOS CCDs using 
an annular region with inner radius 5" and outer radius 
40" to minimize pile-up; the resulting spectra were virtu- 
ally identical to those extracted using the circular region. 
We elected to fit the spectra extracted using circular re- 
gions in order to maximize the photon statistics associated 
with emission features in the Fe K bandpass. 

We inspected the 10-12 keV pn background light curve; 
there were no significant flares (the 10-12 keV rate never 
exceeded 0.04 ct s~^). The total good exposure time was 
90.7 ks for pn and 96.4 ks for each MOS. The 0.2-10 
keV pn light curve, binned to 600 s and normalized by 
its mean count rate (11.74 ct s^^), is displayed in the top 
panel of Figure 1. Similarly displayed in Figure 1 are the 
0.2-1 and 3-10 keV light curves; in these energy ranges, 
the total continuum is dominated by the soft excess and 

The H«/H/3 ratio is intrinsically 3.1, assuming case B recombination and assuming that coUisional excitation is negligible, applicable for 
typical NLR conditions, but see warnings by Netzer 1990 regarding uncertainties in the intrinsic Balmer decrement. 
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hard X-ray power-law components, respectively, as will be 
demonstrated in §4. 

The Reflection Grating Spectrometer (RGS) data were 
extracted using standard extraction procedures, includ- 
ing the SAS task rgsproc and the most recent calibra- 
tion files. We used the first-order data only; spectra were 
binned every 10 channels {'^O.lA, or roughly the RGS res- 
olution). Virtually all bins contained at least 20 counts, 
allowing use of the statistic. Bad RGS channels were 
ignored. The good exposure time was 99.2 ks for each 
RGS. 

The Optical Monitor (OM) was in Science User Defined 
mode, with one "image" window 5.17' x 5.17' centered on 
NGC 3227. The UVWl filter, whose effective area peaks 
near 260 nm, was used throughout. Extraction proceeded 
in a manner similar to §2-3 of Smith & Vaughan (2007). 
Source photons were collected using a 12-pixel radius cir- 
cle centered on NGC 3227. Emission from the host galaxy 
was not subtracted; the effect of including such a constant 
component is that any variability we observe is thus likely 
a lower limit to the intrinsic level of variability associated 
with the AGN. The background was extracted from a cir- 
cular region away from the host galaxy. The resulting net 
light curve, extracted in 1400 s bins, is displayed in Figure 
1. 

2.2. RXTE PCA and HEXTE extraction 

To investigate the form of the >10 keV continuum and 
estimate the strength of any Compton reflection hump 
present, we examined archival Rossi X-ray Timing Ex- 
plorer RXTE Proportional Counter Array (PCA; Swank 
et al. 1998) and High-Energy X-Ray Timing Experiment 
(HEXTE; Rothschild et al. 1998) monitoring data. NGC 
3227 was observed in November 1996 for ~260 ks. NGC 
3227 was also monitored from 1999 January 2 until 2005 
December 4, with one visit every 2-7 days, along with a 
period of more intensive monitoring approximately four 
times daily from 2000 April 2 to 2000 June 7 (see Uttley 
& M'^Hardy 2005 for further details regarding sampling 
patterns). Each monitoring snapshot typically lasted 1 ks. 
There were no observations simultaneous with the 2006 
XMM-Newton observation, so we considered all available 
data in the public archive in order to estimate average 
long-term properties. However, data taken from late 2000 
to early 2001 (modified JuHan day [MJD] 51850 52050), 
affected by the passage of the compact cloud across the 
line of sight, were ignored during spectral fitting. 

PCA data were extracted and screened using standard 
methods and tools. The 'L7-240' background models, ap- 
propriate for faint sources, were used. See e.g., Markowitz, 
Edelson & Vaughan (2003) for details on PCA data ex- 
traction and background subtraction, the dominant source 
of systematic uncertainty (e.g., in total broadband count 
rate) in these data. Counts were extracted from Propor- 
tional Counter Units (PCUs) 0-2, and 2, and 2 only 
for data observed before 1999 December 23, between 1999 
December 23 and 2000 May 12, and after 2000 May 12, re- 
spectively. Response files were generated for gain epochs 
3, 4 and 5 separately (data observed before 1999 March 22, 
between 1999 March 22 and 2000 May 12, and after 2000 
May 12, respectively). As the response of the PCA slowly 
hardens over time due to the gradual leak of propane into 



the xenon layers, data within each gain epoch were further 
split into roughly equal segments, each spanning durations 
of roughly a couple hundreds days, and separate response 
files were generated for each segment. Spectra and re- 
sponses were then added, weighting by exposure times, 
using the FTOOLS SUMPHA and ADDRMF, respectively. 
The total exposure time for all PCA data was 842.7 ks. 

The HEXTE instrument consists of two independent 
clusters (A and B), each containing four scintillation coun- 
ters which share a common 1° FWHM field of view. Source 
and background spectra were extracted from each indi- 
vidual RXTE visit using Science Event data and stan- 
dard extraction procedures. The same good time intervals 
used for the PCA data (e.g., including Earth elevation 
and SAA passage screening) were applied to the HEXTE 
data. To measure real-time background measurements, 
the two HEXTE clusters each undergo two-sided rocking 
to offset positions, in this case, to 1.5° off-source, switch- 
ing every 32 seconds (16 seconds before 1998 Jan.). No 
strong contaminating hard X-ray source within 2° is evi- 
dent (e.g., from the RXTE-VCA or SwiftrBAT all-sky slew 
survey data available at the HEASARC's online SkyView 
service). Cluster A data taken between 2004 Dec 13 and 
2005 Jan 14 were excluded, as the cluster did not rock 
on/off source. Detector 2 aboard Cluster B lost spectral 
capabilities in 1996; these data were excluded from spec- 
tral analysis. Cluster A and B data were extracted sep- 
arately and not combined, as their response matrices dif- 
fer slightly. Deadtime corrections were applied. All data 
within each cluster were combined, except for the 1996 
data; that is, the 1996 (16 s rocking) and 1999-2005 (32 
s rocking) spectra were fit separately. We note that the 
1996 and 1999 2005 spectra agree well with each other 
for each cluster, illustrating HEXTE's performance down 
to source fluxes which are 1% of the background. Good 
exposure times for the four HEXTE source spectra were 
237. f ks (cluster A) and 234.2 ks (cluster B) for the 1999- 
2005 data, and 41.5 (A) and 41.8 ks (B) for the 1996 data. 
All data were grouped as follows: channels 17-30, 31-39, 
40-47, 48-67, 68-79 and > 80 were grouped by factors of 
2,3,4,5,6 and f 0, respectively The LLAGN/LINER source 
NGC 3226 is located about 2' away and is thus in the field 
of view of both the PCA and HEXTE, but its flux is a 
factor of ~ 50 fainter than that of NGC 3227 in the 2-10 
keV band (Binder et al, in prep.; George et al. 2001). 

3. EPIC PN SPECTRAL FITS TO THE FE K BANDPASS 

XSPEC (Arnaud 1996) v.ll.3.2ag was used for all spec- 
tral fitting. All errors below correspond to Ax^ = 2.71 
(90% confidence level for one interesting parameter when 
the errors are symmetric) unless otherwise noted. The 
abundances of Lodders (2003) were used. NGC 3227's 
redshift is 0.00386 (de Vaucouleurs et al. 1991); a distance 
of 20.3 Mpc was used (Mould et al. 2000) was used to 
calculate luminosities. In all models, we included a col- 
umn of neutral absorption fixed at the Galactic column of 
2.1x10^° cm-2 (Dickey & Lockman 1990, Murphy et al. 
1996). 

We started by fitting the EPlC-pn data, 4-10 keV only. 
Residuals to a simple power-law fit {x^ /dof = 1515/999) 
showed an obvious Fe Ka emission line, but also some in- 
teresting structure near 6.0 keV; see Figure 2a. In our fits, 
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we tested two competing models that explain the emission 
features roughly equally well. As explained below, in both 
models, narrow emission features are detected at 6.40 kcV 
and 7.37 keV (Fe Ka and Ni Ka respectively). The mod- 
els differ in how emission near 6.0 keV is modeled: as a 
narrow Gaussian emission line, or as a relativistic diskline 
(hereafter "GA" and "DL" models). 

We started with the former case. First, we included 
a Gaussian component at 6.40 keV to model the Fe Ka 
line. In the best-fit "GA" model, the Fe Ka line best-fit 
energy, intensity and EW were 6.403^Q'Qi[J keV, consis- 
tent with neutral Fe. Its intensity and EW were 3.5±0.4 
xlO~^ ph cm~^ s^^ and 91±10 eV, respectively. Its width 
a was 65±14 eV, which corresponds to a FWHM velocity 
of 7000±1500 km s^^. 

Figure 3 shows contour plots resulting from adding a 
"sliding Gaussian" to the "power-law -I- Fe Ka line" model 
to trace the emission residuals after the narrow Fe Ka line 
was modeled. That is, we added a Gaussian with width a 
fixed at 10 eV and searched over 4-9 keV in units of 0.1 
keV, i.e., in steps smaller than the instrumental resolution. 
The narrow feature near 6.0 keV is clear; other residuals 
near 7.0 and 7.4 keV also warrant further investigation. 

To model Fe K/3 emission, we added a Gaussian to the 
model, keeping the energy centroid fixed and width a tied 
to that of the Fe Ka line. However, the line was not signif- 
icantly detected according to an i^-test. The upper limit 
to the line's EW was 15 eV; the upper limit on the K/3/Ka 
normalization ratio was 0.19 (yielding no useful constraints 
on the ionization state of the Fe-line emitting material). 
We included this feature in the rest of the fits for com- 
pleteness, with intensity fixed to 0.13 times that of the Fe 
Ka line. Residuals to a model consisting of the power-law 
and the Ka and K/3 hnes (where x^/dof was 1028.3/996) 
are plotted in Figure 2b. 

We next added an Fe K edge at 7.11 keV (energy 
fixed) hkely associated with reflection; ^ /dof dropped 
to 1017.6/995; residuals to a model including the edge are 
plotted in Figure 2c. In the best-fit model, the optical 
depth T was 0.05±0.03. 

Next, we used a Gaussian component to model an emis- 
sion hne at 7.39l^;i^? keV, with width a tied to that of the 
Fe Ka line, consistent with emission from Ni Ka. Its inten- 
sity and equivalent width were 41^x10-6 ph cm s and 

13tio eV, respectively, x^/dof dropped to 1010.2/993; it 
was significant at 97.3% confidence according to an F-test 
to include this line. Residuals to a model including this 
line are plotted in Figure 2d. 

Next, we added a narrow Gaussian emission profile near 
6.0 keV to model those residuals; Ax^ was -24.94 for 3 
less do/, significant at >99.99% confidence according a 
standard F-test. This line is not resolved (width a < 200 
eV). Its intensity and EW are 9^4 x 10~^ ph cm s^^ 
and 2lj;g^ eV, respectively. A contour plot of intensity 
versus energy (confidence levels are for two interesting pa- 
rameters) is shown in Figure 4. The best-fit line energy is 
6.041q |!)4 keV inconsistent with 6.24 keV, the lowest pho- 
ton energy associated with a Compton shoulder scattering 
feature. Using a Gaussian to model a Compton shoulder, 

radial emissivity per unit area is quantified as a power law, r~l^ 
13 1 i?g = GMbh/c^ 



we should expect a centroid energy > 6.24 keV; fixing the 
Gaussian line energy at 6.24 keV failed to fully model away 
the 6.0 keV residuals. We performed Monte Carlo simula- 
tions to assess the statistical significance of this emission 
feature, as standard usage of the F-test may overestimate 
the statistical significance of lines at arbitrary energies 
such as this one (see e.g., warnings by Protassov et al. 
2002). The simulations were conducted following §4.3.3 
of Markowitz, Reeves & Braito (2006); see also §3.3 of 
Porquet et al. (2004). The simulations indicated that the 
likelihood that the 6.0 keV emission line is photon noise is 
<0.1%, i.e., the line is significant at >99.9% confidence. 

In this best-fit "GA" model, x^/dof was 983.2/990; F 
was 1.57±0.03. Other best-fit parameters are shown in 
Table 1. Residuals are plotted in Figure 2e. 

Next, we tried to model the 6.0 keV feature as the red 
wing of a relativistic diskline profile. We used a Laor 
model profile (Laor 1991) to model the diskline, with emis- 
sivity index 43 fixed at 3. x^/dof was 988.7/989 (Ax^ = 
-16.8; significant at 99.8% confidence in an F-test) in the 
best-fit "DL" model. The EW of the diskline was Sltfo 
eV ; its inclination was <25°, and the inner radius was <22 
i?J3- Visually, there still remained correlated, emission- 
like residuals from 5.9-6.1 keV, as shown in Figure 2f, 
which might at first suggest that the 6.0 keV emission fea- 
ture may be better modeled as a narrow feature than as a 
red wing of a diskline. However, given the similar values of 
X^/dof and the instrument resolution, it would be difficult 
to demonstrate that the difference in residuals between the 
two models is not consistent with photon noise. 

3.1. MOS 1+2 spectral fits to the Fe K bandpass 

As a double-check on spectral modeling, we applied our 
best-fit "GA" model to both time-averaged MOS spectra 
(fit simultaneously, with all line parameters tied between 
both MOS spectra). The power-law indices and normal- 
izations were allowed to differ between MOS 1 and MOS 
2. Overall, the MOS and PN yielded qualitatively sim- 
ilar results. A simple power-law fit resulted in x^/dof 
= 771.9/558; data/model residuals are plotted in the top 
panel of Figure 5. Adding Gaussian components to repre- 
sent Fe Ka and K/3 emission in a manner similar to the 
pn fits caused x^ jdoj to drop to 546.86/555. Again we 
performed a "sliding narrow [a — 10 eV) Gaussian" test; 
the resulting contours are plotted in Figure 6 and again 
indicate emission near 6.0-6.1 keV. 

Finally, adding a narrow Gaussian (width a fixed at 10 
eV) at 6.1ll[! 'i'^ keV further improved the fit; x^jdof fell 
to 535.87/553 (Ax^ = —11.0 for two less do/; significant 
at 99.7% confidence according to an F-test used in the 
standard way). The best-fit values of the intensity and 
EW were 6.4 ± 3.4 x 10"^ ph cm^^ and 13 ± 7 eV, 
respectively. Data/model residuals for the best-fit model 
are plotted in the bottom panel of Figure 5 (There remain 
~10-15% residuals in ratio space near 7.0 and 7.5 keV in 
Figure 5, but Figure 6 indicates that those features are not 
detected with high significance.). A contour plot of the in- 
tensity of the emission line at 6.11 keV versus line energy 
is plotted in Figure 7. Again we performed Monte Carlo 
simulations to gauge the detection significance, finding the 
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line to be detected at 95% confidence in the MOS 1+2 
spectrum i.e., independent of the pn detection. In fact, 
we can multiply the independent null hypothesis probabil- 
ities for the pn and MOS (< 0.001 and 0.05) to yield a 
"effective" null hypothesis probability of < 5 x 10^'"'. 

The detection of features at 6.0-6.1 keV in simultane- 
ous pn and MOS 1-1-2 spectra illustrates the importance 
of observing with multiple X-ray instruments to attempt 
to distinguish between narrow and/or weak features which 
may be intrinsic to the target and features which arc ar- 
tifacts associated with photon noise (even when Monte 
Carlo simulations are used to gauge the detection signifi- 
cances in individual spectra). The likelihood that a strong 
emission feature could appear at the same energy in both 
the pn and the MOS 1+2 spectra and be due to photon 
noise is likely quite small. 

The photon indices for the MOS 1 and MOS 2 spec- 
tra were 1.37+0.04 and 1.46+0.04, respectively; the flat- 
ter values of the photon index compared to the pn are 
likely the result of pile-up (MOS spectra extracted using 
only pattern events has photon indices which were ~0.3 
steeper than the MOS pattern 0-12 spectra). Finally, we 
added an Fe K edge (energy fixed at 7.11 keV); /dof 
feU to 517.4/552 for t = 0.09 + 0.04. Emission from Ni Ka 
was not significantly detected. All other parameters were 
consistent at the 90% confidence level with those measured 
by the pn. 

Wc also applied the best-fit "DL" model to the MOS 
1+2 spectrum. ^ jdoj was 504.9/550, with all diskline 
and Fe Ka line parameters consistent with those in the pn 
fits. 

4. EPIC PN BROADBAND SPECTRAL MODELING 

Wc started with the best-fit "GA" model and included 
data down to 0.2 keV. In addition to the Galactic column, 
we included a zwabs component at the systemic redshift 
to model any excess cold absorption associated with e.g., 
the host galaxy or circumnuclcar material, A'h, local- 

The residuals (Figure 8a) showed a large soft excess and 
a large absorption trough at ~0.73 to ~0.92 keV which we 
identify as an Fe UTA feature. A narrow absorption fea- 
ture near 1.34 keV is apparent, likely due to Mg XI. A 
narrow absorption feature near 1.85 keV could be due to 
Si XIII, but could also be due to calibration uncertainties 
associated with the instrumental Si K edge. 

We first modeled the soft excess using a steep power- 
law with photon index Fgx near 3. Our final, best-fit 
model (including modeling the warm absorbers and O VII 
emission line; see below) is henceforth referred to as the 
"SXPL (soft X-ray power law)" model. In the best-fit 
SXPL model, Fsx was 3.351q'^q. We also added another 
layer of cold absorption at the systemic redshift, applying 
it only to the hard X-ray power law; in the best-fit SXPL 
model, its column density A^h,hx was 2.9lo'8^10^^ cm~^. 
This improved the overall fit substantially [x^ /dof fell to 
3070.3/1730), but the 0.7-0.9 keV trough remained; see 
Figure 8b. In this model, Fhx was 1.64; we note that 
forcing the photon indices of the soft and hard power-law 
components Fsx and Fjjx to be equal resulted in a very 
poor fit, with x/dof = 6057/1731. 

Next, we used an xstar grid, which assumed a turbu- 
lent velocity width of 100 km s~^, to model a layer of 



warm absorption with ionization parameter log^ and col- 
umn density A^h- The grid assumed an underlying optical 
to X-ray spectral energy distribution described as follows: 
Below 0.001 keV, a power-law component with spectral 
index a = 1.0 (F = 2.0); from 0.001 to 0.04 keV, a power 
law with a = 0.2 (F = 1.2), following e.g., Elvis et al. 
(1994) and Netzer (1996); from 0.04 to 1.0 keV, a steep 
power law with a = 1.9 (F = 2.9; the intrinsic EUV con- 
tinuum is not well studied, but we assume that the soft 
excess seen in the XMM-Newton spectrum extends down 
to 0.04 keV); and above 1.0 keV, a power law with a = 0.5 
(F = 1.5), with the >0.04 keV power law indices based on 
the best-fitting SXPL model. Initially, solar abundances 
and a zero velocity offset relative to systemic were as- 
sumed. Applying one zone of warm absorption to the 
model, x^ldof fell to 2184.1/1728 for log^ ~ 2.2 and ~ 
2x10^^ cm~^. As shown in Figure 8c, however, significant 
data/model residuals remained. A model incorporating 
two zones of warm absorption successfuly modeled all the 
absorption-like residuals, including the trough near 0.75 
keV; x^/dof fell to 1942.1/1726. In this model, the two 
X-ray absorbers had ionization parameters log^io near 0.3 
and log^hi near 2.5. We henceforth refer to the two X-ray 
absorbers as the low-ionization X-ray absorber (though the 
ionization parameter is still higher than that of the opti- 
cal/UV lukewarm absorber; see §9), and high-ionization 
X-ray absorber, respectively. 

By now, the residuals (see Figure 8d) showed a nar- 
row emission feature near 0.56 keV, likely due to O VII. 
Adding a narrow Gaussian component at 0.58+0.01 keV 
with width a fixed at 1 eV caused x^ to drop by 75 for two 
less dof . However, we note that the warm absorbers each 
predict narrow absorption lines due to O VII, so there 
are likely large systematic uncertainties associated with 
the measured intensity of the O VII emission line. This 
model, with xVdo/=1866. 9/1724, is our best-fit SXPL 
model. In this model, the column densities and ioniza- 
tion parameters of the low- and high-ionization absorbers 
are: iVn.io = l-OtJ];? x lO^i cm-^, log^io = lAb+J"^}^, 
A^H.hi = 1.8tJ:^xl02i cm-2, and log^hi = 2.93^°;^^. Other 
best-fit parameters are listed in Table 2. The data/model 
residuals are shown in Figure 8e. Residuals near 1.8 
and 2.2 keV are instrumental (Si K and Au M). Narrow 
emission-like features near 0.42 and 0.92 are tempting to 
identify as N VI and Ne IX emission, but they are likely 
narrower than the instrument resolution. There also ap- 
pear to be some negative residuals near 6.8 keV, close to 
the energies associated with Fe XXV and Fe XXVI. How- 
ever, fitting this feature with an inverted Gaussian does 
not yield significant improvement to the fit; it is likely a 
spurious feature and will not be discussed further. 

Note that residuals near the expected O and Fe L3 edge 
energies due to dust (0.53 and 0.71 keV) are fine; we do 
not require additional edges at those energies (though we 
revisit this issue with the higher resolution RGS data be- 
low). If dust is present in the warm absorber, abundances 
of O and Fe (and also Si and C) may be expected to be 
low due to depletion onto grains (see Snow & Witt 1996). 
However, the measured O and Fe abundances of the warm 
absorber are consistent with solar, and left fixed at solar 
values for the remainder of the paper. 

In the best-fit SXPL model, the hard X-ray power-law 
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component has a photon index Fhx = 1-57 ± 0.02, a cou- 
ple tenths lower than "canonical" values of 1.8-1.9 in many 
broad- line Seyferts. We therefore explored the possibility 
that the observed low value of the photon index is an ar- 
tifact caused by a yct-unmodclcd partial-covering absorb- 
ing component with column density Nn^pc ~ 10^'^ cm~^ 
along the line of sight to the nucleus. Such material is 
expected to exist within light-days of the black hole, as 
demonstrated by the 2000-1 obscuration event. Further- 
more, such a column is required to explain the origin of 
the Fe Ka line (see §9.4) in the absence of strong Compton 
reflection (see §6). Modifying the hard X-ray power-law 
component to be absorbed by a partial-covering neutral 
absorber with iVn^PC set to 0.3, 1.0 or 3.0 xlO^^ cm'^, 
we find only upper limits to the covering fraction in the 
range 5-15%. In the best fit models, Fhx was never higher 
than about 1.6. Partial-covering models are not analyzed 
further. 

As before, we applied our best-fit SXPL model to the 
time-averaged MOS 1-1-2 spectra as a double-check. The 
MOS fits, using this model and fitting over 0.4-10 keV, and 
allowing power-law component parameters to differ be- 
tween MOS 1 and MOS 2, yielded x^/dof = 1406.33/1017, 
and required slightly higher values of the column density 
and ionization parameter of the high-ionization warm ab- 
sorber (TVh.w = l-Stoi X 10^^ cm-2 and logCm = 3.22toTi 
erg cm s^^, respectively) compared to the pn. Fsx for 
MOSl and M0S2 were consistent with that for the pn. 
Fhx for MOSl was 1.52 ± 0.01, slightly flatter than that 
for the pn. Fhx for M0S2 was 1.57±0.02, consistent with 
the pn. All other parameters were consistent at the 90% 
confidence level with those measured by the pn. 

4.1. Alternate parametrizations of the soft excess 

We returned to the pn data to consider alternate meth- 
ods of modeling the soft excess. An alternate parametriza- 
tion of the soft excess, a blackbody, yields a nearly identi- 
cal value of x^/dof (1870.4/1724). Data/model residuals 
are plotted in Figure 8f. The blackbody temperature ksT 
= 83tl cV in the pn fit (77 ±3 eV for the MOS), similar to 
values found for model fits incorporating blackbody com- 
ponents for most other Seyferts, e.g., ksT ^100-150 eV 
for many narrow-line Seyfert Is. It is generally accepted 
now that a blackbody is a physically unplausible descrip- 
tion of the soft excess in AGN. Bechtold et al. (1987) were 
the first to point out that the best-fit blackbody tempera- 
tures are generally too high to be associated with accretion 
disks around supermassive black holes. The consistency 
of blackbody temperatures across a wide range of Seyfert 
1 properties, including black hole mass, furthers suggests 
that a blackbody is an unphysical parametrization (e.g., 
Gierlihski & Done 2004). However, the best-fit blackbody 
temperature in the current fit is consistent with previous 
results for NGC 3227 (e.g., Komossa & Fink 1997b). We 
note that in this model (henceforth denoted simply as the 
"BB" model), A^h.hx bas dropped to < A^h, local, i-e., the 
value of A^H,HX is model dependent. 

We next tried fitting the soft excess using the Comp- 
tonization model CompST (Sunyaev & Titarchuk 1980). 
Such components have been used e.g., by Gierlinski 
& Done (2004) to model the soft excess emission as 
Comptonization of accretion disk seed photons in a cool 



(k-BT ^ 0.3 keV), optically-thick corona, distinct from 
the optically-thin, hot (fceT ~ 100 keV) corona gener- 
ally thought to be responsible for the hard X-ray power- 
law component in Seyferts. Physically, such a component 
could be potentially associated with the transition region 
between an optically-thick accretion disc and an optically- 
thin inner corona (Magdziarz et al. 1998), or in the hot, 
ionized surface of the accretion disk (Hubeny et al. 2001; 
Janiuk, Czerny & Madejski 2001). However, similar to the 
blac;kbody component, using such Comptonization c;ompo- 
nents to model soft excesses tends to yield a narrow range 
of best-fit temperatures across a range of Seyfert proper- 
ties. In our best-fit model ("COMPST", where x^/do/ 
= 1863.4/1723), ksT = 0.35to°ol keV, and optical depth 
T = 24^4, similar to previous fits to other AGN (however, 
see e.g., Vaughan et al. 2002 for warnings regarding the co- 
variance of optical depth and temperature in this model). 
Best-fit parameters for the models using the blackbody 
and CompST components are listed in Table 2, though 
the reader must bear in mind the difficulties each of these 
models faces in terms of physical plausibility. Data/model 
residuals are not plotted, as they are virtually identical to 
those for the SXPL model (Figure 8e). 

Figure 9 shows unfolded model spectra for the SXPL, 
BB and COMPST models, demonstrating how the steep 
soft excess is produced in each (and demonstrating why the 
value of A^H.HX depends on how the soft excess is modeled). 

Many recent forays into modeling soft excesses have 
found success using ionized reflection models modified by 
relativistic smearing, i.e., a blurred ionized reflector (e.g.. 
Crummy et al. 2006); an origin in atomic features could 
plausibly explain the consistency of soft excess features 
across a range of Seyfert black hole masses. Starting with 
the best-fit SXPL model, we removed the soft X-ray power- 
law and the 6.04 keV emission line, and added a Ross 
& Fabian (2005) reflection component convolved with the 
kernel associated with relativistic motion around a Kerr 
black hole (a Laor profile). The photon index of the il- 
luminating continuum was tied to that of the hard X-ray 
power-law. The disk outer radius was fixed at 400 Rg] 
the emissivity index was fixed at 3. The 6.4 keV Fe Ka 
and Ni Ka line parameters were fixed at the values in the 
time-averaged spectrum. The best-fit model had x'^/dof 
= 2698.4/1721, with large (±10%) correlated residuals be- 
low ~2 keV (see Figure 8g) . The best- fit Fe abundance was 
<0.15, and the best-fit value of the reflector ionization pa- 
rameter ^ pegged at the lower limit of 30 erg cm s~^, likely 
indicating that lower values of ^ may be more appropri- 
ate. The low level of Compton reflection > 10 keV in NGC 
3227, R < 0.5, as demonstrated in §6, also suggests that 
ionized reflection in NGC 3227, if it exists, is not strong. 
We do not consider this model further. 

Finally, we tried to model the soft excess using smeared 
absorption to model ionized absorbing material moving at 
relativistic velocity, such as a wind launched from the in- 
ner accretion disk SWiNDl (Gierlinski & Done 2006). As- 
suming a single power-law to model the intrinsic contin- 
uum emission, modifled by smeared absorption, yielded 
a poor fit {x /dof = 8388/1728). Modeling the intrinsic 
continuum using two power-laws, as in the SXPL model, 
and adding smeared absorption plus the low-ionization X- 
ray absorber, yielded a good fit {x^/dof = 1892/1723), 
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in which a smeared ionized absorber with log^ ~2.7 erg 
cm s~^ was able to mimic the high-ionization X-ray ab- 
sorber, given the resolution of the pn. However, analysis 
of the RGS spectrum (§5) supports the need for both the 
low- and high-ionization X-ray absorbers, and so we do 
not consider this model further. 

5. RGS MODELING OF THE IONIZED ABSORBERS AND 

EMISSION LINES 

We first modeled the continuum of the RGS spectrum 
with a model of the same form as the SXPL model from 
the pn fits: a steep soft power-law with Fsx ~ 3.2 absorbed 
by a cold column A^h, local = 1 x 10^^ cm~^, and a hard 
power-law absorbed by a column A^h.hx ^ 1 x 10^^ em~^, 
and dominating the continuum only above ~ 1.5 keV (and 
with both A^H.HX and Fhx very poorly constrained, given 
that the RGS only covers up to 2 keV). The residuals to 
this model {x^ /dof = 1156/433) are shown in Figure 10. 

There is a broad absorption feature at ^740 780 eV, 
which we identify as an Fe M-shell UTA trough. Its en- 
ergy range (assuming a velocity offset from systemic near 
zero) indicates that log^ of the gas where this feature orig- 
inates must be in the range ~ -1-0.3 to -1-1.4, according to 
an XSTAR table; and absorption is mainly due to species 
of Fe in the ionization range ^ VI- XIII, and likely not 
due to Fe ~ XIV- XVI (Gu et al. 2006) 

Many narrow absorption features are evident; the ten- 
tative identifications (and lab-frame energies) of the most 
visually prominent ones are labeled in Figure 10 and in- 
clude C VI Lya (368 eV), C VI hy(3 (436 eV), N VII (500 
eV), O VIII Lya (653 eV), O VII He/3 (666 eV; perhaps 
blended with the N VII edge), O VII He7 (698 eV), Fe XVII 
3d-2p (812 eV), Ne IX (922 eV, though this line may be 
blended with Fe absorption lines at ~ 918-934 eV), Ne X 
(1022 eV) and Mg XI (1352 eV). 

A strong edge at 739 eV due to O VII is not obvious, 
though it could be blended with the Fe UTA feature. Sim- 
ilarly, an edge at 870 eV due to O VIII is not obvious. The 
lack of a very strong (> 20 — 30% drop) O VIII edge sug- 
gests that any absorbing material containing O VIII has a 
hydrogen column density TVh < 10^^ ''' or so. There are also 
narrow emission features near the energies for the O VII 
and N VI resonance lines (574 eV and 428 eV, lab- frame); 
these are discussed further below. 

We next modeled the two ionized X-ray absorbers dis- 
cussed in the pn fits above, using the same xstar ta- 
ble models, with their respective velocity offsets relative 
to systemic initially frozen at zero, and with abundances 
fixed to solar values, jdoj fell to 955.40/429 for val- 
ues of log^io and log^hi near 1.2 and 2.8, respectively, and 
with column densities for each near 1 x 10^^ cm~^. How- 
ever, many of the absorption lines seemed to be blueshifted 
slightly. 

We tried various combinations of allowing the absolute 
redshift of each absorber, z\o and Zhi, to be thawed from 
the sytemic value. Thawing 2:hi only yielded y^/dof = 
884.87/428 for a best-fit value of zu near -0.0036 (-0.0075 
relative to systemic). Thawing both Zhi and z\o, but keep- 
ing their values tied, yielded y^/dof = 887.46/428 for a 
best-fit value of Zhi = zio = -0.0033 (-0.0072 relative to 
systemic). Finally, thawing both Zhi and z\o and allowing 
them to vary independently, yielded /dof = 864.09/427 



for best-fit values of zm — -0.0034 (-0.0073 relative to sys- 
temic) and z\o = -1-0.0016 (-0.0023 relative to systemic). 

Finally, to this last model (which had the lowest 
value), we added five narrow (width cr=0.5 eV) Gaussian 
components to represent emission at the systemic redshift 
from the O VII (/), (i) and (r), N VII, and N VI (r) 
(note that we model the sum of the two intercombina- 
tion lines for O VII). Their measured energies, intensi- 
ties, and fit parameters are listed in Table 3. N VI (/) 
and N VI (i) emission lines were detected only as upper 
limits, with intensities < 6.7 x 10~^ and < 3.0 x 10"^ 
ph cm~^ s~^, respectively, and so are not listed in Ta- 
ble 3. However, as for the pn fits, we caution that since 
both warm absorbers require some O VII, N VI and N VII 
absorption, the emission line intensities likely have large 
systematic uncertainties. Final parameters for the ion- 
ized absorbers in our best-fit SXPL model (with x^/dof 
= 785.7/417) are log6o = l-2lt°om, Nu,io = x lO^i 

cm-2, logai = 2.90l°:2i, and A^h.m = 2Atli x lO^i cm'^; 
other model parameters are listed in Table 4. We note that 
Fsx = 3.00 ± 0.25 and iVH.iocai = lO.S+ifj x 10^" cm-^ 
both consistent with the best-fit SXPL model for the EPIC 
data. 

In this model, the best-fit value for the red- 
shift of the high-ionization absorber is Zhi (absolute) 
= -0m302tomall or Zhi (relative to systemic) = 
— 0.00688^Q Q0080' which corresponds to an outflow velocity 
of -(2060^^7q) km s~^ relative to systemic. For the low- 
ionization absorber, zio (absolute) = -1-0.002461" ;QQQg4 or 
zio (relative to systemic) = — 0.001401q;qoo64! correspond- 
ing to an outflow velocity of -(420l^9o) ^™ relative 
to systemic. Given the best-fit values of the two warm 
absorbers here and the relative line strengths expected in 
those cases, the outflow velocity for the high-ionization 
absorber is constrained primarily by the narrow absorp- 
tion lines due to Ne IX, Ne X, and Mg XI; constraining 
the respective energies of each of these lines separately 
yielded corresponding outflow velocities relative to sys- 
temic of 1300t^^2, 600 ± 600, and 1800 ± 600 km s-\ 
respectively. Isolating the observed C VI Ly/3, N VII and 
O VIII Lya absorption lines yielded similar velocity off- 
sets (2700 ±700, 1800 ±600 and 1840 ±460, respectively), 
though contributions to these lines from both the low- 
and high-ionization absorber are expected, given the best- 
fit ionization levels. The only strong absorption features 
which are unique to only the low-ionization absorber and 
unambiguously studied here are the (possibly blended) Fe 
UTA feature and the O VII edge. Models wherein the out- 
fiow velocity of the low-ionization absorber is either much 
closer to systemic, equal to that of the dusty lukewarm UV 
absorber (tens of km s^^; Crenshaw et al. 2001), or equal 
to that of the high-ionization absorber are not definitively 
ruled out at high confidence. Setting the absolute value 
of zio equal to systemic, +0.00353 (-0.00033 relative to 
systemic, corresponding to an outfiow velocity of 100 km 
s~^), or Zhi or +0.003859 and refitting yielded increases 
in of only 9, 6, and 16 respectively. In each case, the 
outfiow velocity for the high-ionization absorber remained 
virtually unchanged from the best-fitting SXPL model. 

As an additional confirmation on the xstar table re- 
sults, we measured the EWs of selected individual ab- 
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sorption features. The measured EW for the Fe UTA (- 
(ISj^g)) eV is consistent with the prediction from an XSTAR 
tabic with A^h.Io = 1-1 x lO^^ cm'^ and log6o=l-21. The 
measured EWs for Nc IX (-5.5±2.0 cV), Ne X (-4.4±1.6 
eV) and Mg XI (— 4.1 ± 2.0 eV) are in good or reason- 
able agreement (given the data quality) with the pre- 
dicted EWs from an XSTAR table with iVH,hi = 2.4 x 10^^ 
cm~^ and log^hi=2.90 (these lines are not expected to be 
generated with significantly noticeable EW in the low- 
ionization absorber). The measured EWs for C VI Ly 
13 (-1.0 ± 0.9 cV), N VII (-(l.lto:2) eV) and O VIII 
(— 3.3±0.7 cV) are in good or reasonable agreement (given 
the data quality) with the sums of the predicted EWs from 
the two XSTAR tables. 

Most of the remaining data/model residuals are con- 
sistent with RGS calibration uncertainties or arc instru- 
mental in nature (Pollack et al. 2007; observed energies 
listed): a small dip near 394 eV; large residuals near the 
O instrumental edge, 524-530 eV; a small dip at 573-577 
eV (which further adds to the systematic uncertainties as- 
sociated with O VII emission lines intensity); residuals at 
~680-695 eV, associated with an instrumental fluorine fea- 
ture; an apparent edge at 850 eV; and a dip at 937-943 
eV. Many of these instrumental features are also seen in 
other RGS spectra of Scyfcrts (e.g., Braito et al. 2007). 

Wc rely on the RGS primarily for determination of warm 
absorber and soft X-ray emission line parameters; all hard 
X-ray parameters are best constrained by the EPIC data. 
Still, for the SXPL model, most best-fit parameters (in- 
cluding those for the warm absorbers) were consistent at 
the Ax^ = 2.71 significance level between the RGS and 
the pn; A?h,hx and Fhx were the exceptions, and clearly 
better constrained by the EPIC. Simultaneous RGS/pn fit- 
ting using the SXPL model {x^/dof = 3220/2161) yielded 
best-fit values for Fsx, .^h, local and warm absorber column 
densities and ionization parameters consistent with those 
measured by the RGS and EPIC-pn separately. Similarly, 
modeling the continuum in the RGS data using a black- 
body or a CompST component yielded good fits {y^ /dof 
= 805.7/417 and 806.4/416, respectively), with a best-fit 
blackbody temperature and CompST parameters consis- 
tent with those measured using only EPIC data, and with 
no significant changes to the warm absorber or emission 
line parameters compared to using a soft power-law com- 
ponent to model the soft excess. 

To this point, we have been assuming that the ionized 
absorbers along the line of sight have been full-covering. 
Using the RGS data, we tested a model wherein only a 
fraction of the continuum radiation is absorbed by the ion- 
ized absorbers. We found the upper limit to the fraction of 
continuum radiation that remains unattenuated by ionized 
absorption to be 30%, i.e., a covering fraction of > 70% 
for the RGS data. For the pn data, the covering fraction 
is inferred to be > 93%. 

We also tested for the presence of radiative recombina- 
tion continuum (RRC) features in the RGS data associated 
with H- or He-like ions, expected if the emission lines are 
associated with gas that is photo-ionized. Using a REDGE 
component in XSPEC with energy fixed and with k^T fixed 
at 0.5 keV (arbitrary value), we found upper limits in the 
range 5-8 eV for C V, C VI, N VI and O VII RRCs, upper 
limits < 20 eV for O VIII and Ne IX, and an upper limit 



near 50 eV for an Ne X RRC. A N VII RRC at 667 eV 
significantly improved the fit, but that was due to fitting 
the instrumental fluorine feature, and therefore is likely 
not real. Finally, we note no significant emission from the 
Fe XVII L line (3d 2p) at 826 cV (upper limit of 1 eV), 
limiting the possibility that the observed emission lines are 
generated via coUisional ionization processes. 

We now discuss the search for X-ray spectral features 
in the RGS associated with dust embedded in the warm 
absorbers. The expected effects of embedded dust on the 
X-ray spectrum are: 1) Edges due to neutral Fe, O, C, Si 
and Mg (depending on dust composition) in excess from 
the absorption expected from gas. The O edge would be at 
531 eV. The three Fe edges are the L3, L2, and LI edges, 
at 707, 721 and 846 eV (assuming oxidized Fe; edges due 
to pure Fe would be ~3 eV higher), roughly in a 10:5:1 
ratio (e.g., Kortright & Kim 2000; Bearden & Burr 1967; 
Schulz et al. 2002). 2) There may be weak, narrow absorp- 
tion lines due to Fe oxide species, in particular at 702 eV 
(e.g., Crocombette et al. 1995). 3) Fe, C and O abundances 
in the warm absorber would be low due to depletion onto 
dust grains (Snow & Witt 1996; Komossa & Fink 1997a). 

Near 707 cV and 721 cV (the Fe L3 and Fe L2 edge 
energies) the residuals in the RGS spectrum of NGC 3227 
do not appear obviously edge-like. Inserting an edge at 
any of neutral Fe or O energies does not improve the fit; 
upper limits to the optical depths of O, Fe L3 and Fe L2 
in excess of the neutral absorption associated with the gas 
components modeled above (both local and Galactic) are 
0.17, 0.05 and 0.21, respectively. In contrast, the Fe L3 
edges seen in Cyg X-1 (Schulz et al. 2002) and MCG- 
6-30-15 (Lec ct al. 2001) with the Chandra HETGS are 
both >50% drop across the edge (r > 0.7). From the Fe 
L3 edge in NGC 3227, the implied upper limit on the Fe I 
column density is ~ 10^'' cm~^. Assuming the abundances 
of Lodders (2003), the corresponding equivalent hydrogen 
column density is constrained to be less than ~ 10^^'^ 
cm~^. This limit is consistent with the equivalent hydro- 
gen column density implied by optical line absorption (see 
§1). The limit on the EW (relative to locally-determined 
and unabsorbed continuum) of a narrow (width a ~ 0.5 
eV) absorption line at 702 eV is -0.5 eV. Finally, thawing 
the abundances of Fe, O, and Mg in the warm absorbers 
(abundances in both absorbers tied) does not yield signif- 
icant evidence for deviations from solar abundances. 

6. CONSTRAINING THE >10 KEV CONTINUUM WITH 
RXTE AND SWIFT-BAT 

We fit archival RXTE data and Swift Burst Alert Tele- 
scope (BAT) 9-month survey data to constrain the amount 
of the Compton reflection and constrain any high-energy 
power-law cutoff. We did not do simultaneous fitting with 
the EPIC spectrum, since the XMM-Newton and latest 
RXTE and BAT data obtained were separated by a year 
and since the only two broadband model componcnits in 
the range of overlap are the variable hard X-ray power-law 
and the hard X-ray absorber, whose column A^h,hx ^ 10^^ 
cm~^ is not well constrained by the PGA. However, the 
observed 2-10 keV flux during the 2006 XMM-Newton ob- 
servation, i^2-io = 3.5x10"^^ erg cm~^ s~^, fell in the 
approximate range probed by the RXTE monitoring; only 
during MJD '~51850-52050 was i^2-io observed to be con- 
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sistently below 2x10^^"'^ erg cm"^ s~^. This suggests that 
XMM-Newton caught the source in a "typical" flux state. 

We fit the PCA data over the 3.3-30 keV range; the 
source is not well detected at higher energies. HEXTE 
data below 24 keV were ignored as the responses of clus- 
ters A and B diverge somewhat for faint sources (N. Sha- 
poshnikov, private communication, 2005); there is good 
agreement above this energy. HEXTE data above 100 
keV were also ignored. 14-195 keV BAT data were in- 
cluded; four-channel spectra and response matrices were 
taken from the BAT 9-month AGN catalog websitcEl- A 
constant coefficient was included in all spectral models to 
account for minor normalization offsets between the PCA 
and HEXTE; its value for the HEXTE spectra relative to 
the PCA spectrum was typically 0.8-0.9. No constant was 
needed for the BAT spectrum relative to the PCA spec- 
trum. 

We applied a model consisting of a power-law (initially 
with a high-energy cutoff set arbitrarily at Ec 500 keV), 
a Gaussian emission line at 6 keV to model Fe Ka emis- 
sion, and a component to model Compton reflection from 
neutral material, using pexrav (Magdziarz & Zdziarski 
1995), all modified by a cold column A'h.hx constrained 
to be less than 6 x 10^^ cm~^. The pexrav component's 
inclination was fixed at 45°, solar Fe abundances were as- 
sumed, and the high energy cutoff was tied to that of the 
power-law. The best- fit model had /dof — 297/111, the 
photon index Thx was 1.66±0.01, and the value of the 
Compton reflection strength was 0.16^g'Q2. Allowing 
the high-energy cutoff to be free and rc-fitting, /dof fell 
to 222/110 (significant at 5.6cr confidence according to an 
F-test). Thx was 1.63±0.02, R was 0.40±0.07, and ^ 
90±20 keV. The cold column Nh,hx was G.OljJ § x lO^i 
cm~^ (uncertainty pegged at upper limit). The residu- 
als to the best-fit models with and without a high-energy 
cutoff are shown in Figure llF^. 

Because the RXTE and XMM-Newton observations 
were not simultaneous, it is not clear whether the slightly 
flatter value of Fhx for the pn in the SXPL model 
(1.57±0.02) compared to that for the time-averaged RXTE 
spectrum is due to source variability or associated with 
systematic differences in modeling spectra taken with two 
instruments with very different energy resolution. For ex- 
ample, the average Fe Ka line flux as measured by the 
PCA was 9 ± 1 X 10~^ ph cm~^ s~^, but this is not au- 
tomatically an indication of line variability between the 
RXTE and XMM-Newton observations given the differ- 
ences in PCA and EPIC resolution, the systematic uncer- 
tainties associated with fitting such a complex model to 
the low-resolution PCA data, and possible uncertainties 
associated with PCA/pn cross-calibration. All uncertain- 
ties listed here are one-dimensional statistical errors only 
and do not account for systematic uncertainties associ- 
ated with correlations between parameters or for system- 
atic instrumental uncertainties associated with observing 
faint objects. However, we can safely conclude that the 

http:/ /swift. gsfc.nasa.gov/docs/swift/results/bs9mon/ 

R = f2/27r, where Vt is the solid angle subtended by the reflector 

The presence of the high-energy cutoff is also inferred from fitting PCA+HEXTE only: excluding the BAT data, and refitting, we find that 
including a high-energy cutoff still significantly improves the fit, with best-fit values consistent with the PCA -|- HEXTE + BAT fit, though 
the statistical uncertainty on Ec is 1.8 times larger. Ignoring the HEXTE data, and fitting PCA + BAT only, the best-fit values of Vjjx, R, 
and Ec are 1.62ljj °4, 0.53^° and esljg keV, respectively. 



value of R must be relatively low: Figure 12 shows con- 
tour plots for Ec and i? as a function of Thx from the 
joint RXTE-\-Swift-}i AT fit, demonstrating that values of 
i? > 0.6 are rejected at high confidence. 

7. TIME- AND FLUX-RESOLVED SPECTRAL FITS 

We investigated the temporal behavior of the X-ray 
spectrum of NGC 3227 on a range of time scales. Any 
rapid variation in ionization parameter would rule out 
models where the warm absorber is spatially extended 
and/or located very far from the black hole; tracking ion- 
ization responses to continuum flux variations can yield 
constraints on the location of the warm absorbers (e.g., 
Krongold et al. 2007) . Similarly, any rapid variation in the 
strength of an Fe K emission feature (e.g., Tombesi et al. 

2007) would suggest an origin close to the X-ray illuminat- 
ing source. We used the XMM-Newton spectrum to study 
variability on time scales < 1 day (hereafter "short" time 
scales); we used the i?A'T'i? intensive monitoring from 2000 
April to June to investigate variability on time scales from 
^2 weeks to 2 months (hereafter "medium" time scales), 
and the 1999-2005 monitoring to investigate variability 
on time scales from months to several years ("long" time 
scales"). 

7.1. XMM-Newton time-resolved spectral fits 

We investigated the < 1-day time-resolved spectral be- 
havior of NGC 3227 during the XMM-Newton observation 
by first splitting the EPIC-pn data into five segments of 
equal duration (good exposure time 18 ks); shorter dura- 
tions would not have yielded adequate photon statistics in 
narrow features. In the Fe K bandpass, spectral fitting re- 
vealed several "candidate" narrow transient emission fea- 
tures between 4.5 and 5.2 keV, each with Ax^ in the range 
-11 to -14. Monte Carlo simulations (see §3) showed that 
each of these emission features were significant at 92-97% 
confidence. However, when one takes into consideration 
the number of segments searched over (Vaughan & Uttley 

2008) , the significance levels drop to 60-85% confidence; 
these features are not discussed further. Considering the 
broadband spectrum, all parameters for the neutral and 
ionized absorbers and the O VII emission line were consis- 
tent with those found in the time-averaged spectrum. Sim- 
ilarly, flux-resolved spectral analysis, performed on high- 
and low-flux pn spectra (extracted above and below a 0.2- 
12 keV pn count rate of 11.0 ct s~^, respectively), yielded 
no significant evidence for variability of the neutral or ion- 
ized absorbers, the 6.4 keV Fe Ka line, or the O VII line. 
The 6.04 keV line was not significantly detected with these 
short exposures. 

However, we noticed in these time-resolved fits that the 
normalization of the soft excess increased by ~25% from 
the first 20 ks segment to the second segment. Similar 
results were obtained when modeling the soft excess as a 
steep power-law, a blackbody, or a CompST component. 

The rapid (ks and shorter) variations seen simultane- 
ously in both the hard and soft X-ray light curves (Figure 
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1) are similar to those routinely detected in other Seyferts; 
the associated variability mechanisms will be explored in 
depth in a future paper (Arevalo et al., in prep) and are 
not discussed further here. Here, we are concentrating on 
the soft X-ray band's relatively slower and distinct trend. 
Since this trend occurs over tens of ks, we will refer to it 
as "rapid" for the remainder of this paper. 

We repeated the time-resolved analysis, dividing the 
data into ten segments of equal duration (9 ks good expo- 
sure time), fitting the SXPL model, and keeping all warm 
and cold absorber, O VII line and Fe emission line param- 
eters frozen at their time-averaged values. 

The resulting light curves of Fhx, Fgx and the normal- 
izations of the SXPL and HXPL components are plotted 
in Figure 13. The SXPL normalization seems to track the 
light curve of 0.2-1 keV flux in Figure 1. To further illus- 
trate the rapid variability of the soft excess, the data for 
segments 1 (the lowest soft excess flux bin), 2, and 6 (the 
highest soft excess flux bin) are plotted in Figure 9. 

Finally, fractional variability amplitudes (i^var see 
Vaughan et al. 2003 for a definition) were calculated for 
16 energy bands. The resulting i^var for the entire dura- 
tion, for the first 20 ks, and for the final 80 ks are plotted 
in Figure 14, further illustrating the rapid variability in 
the soft excess. Further detailed variability analysis of the 
variable X-ray continuum, including coherence, intra- X- 
ray time lags and power spectral density function mea- 
surement, will be presented in a future paper (Arevalo et 
al., in prep). 

7.2. RXTE time-resolved spectral fits 

Time-resolved spectral fitting of the RXTE data closely 
followed Markowitz, Edelson & Vaughan (2003). Bin sizes 
were chosen to optimize the trade off between minimiz- 
ing the uncertainties on Fe line flux and maximizing the 
number of bins. This yielded 5 bins, each of duration 13 
days, for the medium time scale. On the long time scale, 
start /stop times were chosen to avoid 60-days gaps where 
RXTE did not observe the source due to sun-angle con- 
straints, with most durations roughly 100 days before loss 
of PCUO on 2000 May 12, and 150 days afterward, yielded 
16 bins. Response matrices were generated separately for 
each segment. 

We used a model consisting of a power law component, 
a Gaussian component to model Fe K emission, and a 
PEXRAV component to model Compton reflection, all ab- 
sorbed by cold material (with a column density A'^h.hx con- 
strained to be < 6 X 10^^ cm^^, unless the bins included 
data taken during the 2000-1 obscuring event). In the best 
fit models, the best-fit values of R, Fe line rest-frame en- 
ergy centroid, and Fe line width a were OAtoi (0.2±0.1), 

6.261°;}? (6.27l°:°D keV, and O-bOt^il (0.331°;°^) keV 
in the medium (long) time scale fits, respectively. To de- 
termine which parameters to leave free during each time- 
resolved fit, we first fit all segments simultaneously with all 
parameters tied, and thawed one parameter at a time, test- 
ing for significant improvement in the fit according to an 
F-test. We found it was significant to thaw the power law 
normalization, Fhx, Fe line intensity /fcKq, and A'h.hx in 
the individual, time-resolved fits. 

1^ See §3.3 of Markowitz, Edelson & Vaughan (2003) for further det 
calculated at http://www.wessa.net/rankcorr.wasp 



Errors for /fcKq and Fhx were derived using the point- 
to-point varianceP^. Errors on values of ^2-10 within a 
time bin were determined from the 2-10 keV continuum 
light curve, using the mean flux error on the ~1 ks expo- 
sures in that time bin. Figure 15 shows the resulting light 
curves for ^2-10, hcKa, and Fhx- 

Fractional variability amplitudes -FVar were calculated. 
For the medium time scale, Fvm- for the i^2-io and /peKa 
hght curves were 38.4±0.3% and 22.7±10.3%, respectively, 
although we caution that with only 5 data points, these 
measured values of -FVar are likely not highly reliable. 
For the long time scale, Fvar for F2-10 and IpeKa were 
32.7±0.7% and 20.0±4.3%. Qualitatively, these results 
are consistent with what Markowitz, Edelson & Vaughan 
(2003) found for a small sample of Seyferts: the Fe line 
does not vary as strongly as the continuum flux. 

Inspection of the F2-10 and /fcKq light curves would 
seem to indicate, to the human eye at least, similar vari- 
ability trends. However, with so few points (especially for 
the medium time scale), such a statement is not highly 
significant. On the medium time scale, the best-fit lin- 
ear relation is /fcKq = 2.05 ± 1.31 x ^2-10 + 4.1 ± 3.1, 
with IpeKa in units of 10~^ ph cm~^ s~^ and F2_io in 
units of 10~^^ erg cm~^ s""'^. Figure 16 shows the zero-lag 
correlation plot of iFeKa as a function of i^2-i0i with the 
best-fit linear relations plotted as dashed lines. The zero- 
lag Spearman rank correlation coefficient is 0.700; the 
corresponding null hypothesis probability (the probabil- 
ity that the correlation could arise from randomly chosen 
data points) is 16%; i.e., the correlation is significant only 
at 84% confidence. 

For the long time scale data, the best-fit linear relation 
is hcKa = 1-44 ± 0.30 X i^2-io + 3.5 ± 1.2. The zero-lag 
Spearman rank correlation coefficient is r=0.739, signifi- 
cant at 99.6% confidence. We searched for lags using an 
Interpolated Cross Correlation Function (ICE; Gaskell & 
Peterson 1987, White & Peterson 1994), with errors de- 
termined using the bootstrap method of Peterson et al. 
(1998). The ICE peak correlation coefficient rmaa;=0.854 
was reached at a delay of 75 ± 690 days {F2-10 leading 
IpcKa), i-e., consistent with zero lag. 

With such few data points and low correlation coeffi- 
cients on both medium and long time scales, any claim of 
a correlation must be deemed tentative at best. Cross cor- 
relation analysis is further complicated by the red-noise 
nature of the continuum (and likely line) light curves. 
Specifically, cross-correlation between two unrelated red- 
noise light curves can randomly yield spurious correla- 
tions with higher than expected values of the correlation 
strength r^ax (e.g., Welsh 1999). Lags where rmax is not 
very close to 1.0 should thus be treated with skepticism. 
Additional monitoring, spanning much longer durations 
and wider ranges in both -F2-10 and /fcKq and encom- 
passing additional upward/downward trends in the light 
curves of both F2-10 and /fcKo are required to critically 
test for any significant correlation. We therefore conclude 
that there is, at best, tentative evidence from the RXTE 
time-resolved spectral fits for a significant fraction of the 
the Fe line intensity to respond to continuum variations 
on time scales shorter than 700 days. The lack of strong 

and the definition of the point-to-point variance 
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hard X-ray variability during the XMM-Newton observa- 
tion (i^var = 8.5±0.2% for the 3-10 keV Ught curve, binned 
to 600 s) means we cannot draw any conclusions about the 
response of the line on < 1 day time scales. 

8. x-ray/uv continuum light curve correlations 

The OM light curve of NGC 3227 displayed in Figure 
1 shows a ~10% increase across the observation, with an 
RMS of 2.8%. Rapid optical/UV variability in Seyferts 
on time scales of < 1 day is relatively rare. Of the sam- 
ple of 8 Seyferts examined by Smith & Vaughan (2007), 
NGC 3783 displayed the strongest variability on < 1 day 
time scales, with an RMS of 2.9% in the UVW2 filter; 
there were only 3 additional observations using either the 
UVW2 or U filter where the RMS was > 1.8%. The cur- 
rent XMM-Newton data would thus seem to indicate that 
NGC 3227 is towards the top of the list of Seyferts which 
display strong UV continuum variability on < 1 day time 
scales. In NGC 3227, the X-rays are much more strongly 
variable on these time scales (RMS for the 0.2-1 and 3- 
10 keV light curves were 15.8% and 8.9%, respectively), a 
result similar to what Smith & Vaughan (2007) found for 
their sample. 

Visually, it is tempting to connect the gradual bright- 
ening observed in the OM light curve with that observed 
simultaneously in the soft X-ray band during the XMM- 
Newton observation. However, the rise in soft X-ray flux 
(~40%) is much greater than that observed in the UV 
band (^10%). We calculated ICFs and Discrete Corre- 
lation Functions (DCF; Edelson & Krolik 1988) between 
the UV and soft X-ray light curves (the latter rebinncd to 
1400 s); the results are plotted in Figure 17, along with 90 
and 95% confidence limits from the Bartlett method. We 
found no significant correlation: rmax peaks at only ^0.5; 
there are no "bends" or multiple trends in the OM light 
curve to drive a correlation. The warning regarding cross- 
correlation analysis performed on two unrelated red-noise 
light curves bear repeating: beware of spuriously high val- 
ues of Tniax (e.g., Wclsh 1999) and treat values of r^ax 
which are not very close to 1.0 with skepticism. The sug- 
gestion of a correlation between the soft X-ray band and 
the UV band, while visually tempting, is thus speculative 
at best. 

9. DISCUSSION 

We have analyzed a ~100 ks XMM-Newton long-look 
of NGC 3227, observed in December 2006. In both the 
EPIC-pn and RGS spectra, we have modeled the ionized 
X-ray absorber using two components, with very similar 
column densities, and with ionization parameters log^ near 
1.2 and 2.9 (in the RGS spectrum). With the RGS, we 
have constrained the outfiow velocity of the high-ionization 
absorber to be -(2060ti7o) ^™ relative to systemic. 
The best estimate for the outflow velocity relative to sys- 
temic for the low-ionization absorber is (420l^gQ) km s~^, 
though a wider range of possible outflow velocities cannot 
be ruled out at very high confidence. In §9.1, we discuss 
further details of these ionized outflows, including explor- 
ing connections to the low-velocity, very low ionization, 
dusty lukewarm absorber observed in the UV. 

The steep soft excess, which does not seem be to affected 
by the same neutral material obscuring the hard X-ray 



continuum, is shown with the EPIC pn to be rapidly vari- 
able compared to other Seyferts on < 1 day time scales; its 
normalization increases by ~25% in ~20 ks. The UV con- 
tinuum light curve during the XMM-Newton observation 
also shows a relatively strong increase compared to most 
other Seyferts on < 1 day time scales. Possible origins for 
these behaviors are discussed in §9.2. 

The hard X-ray power-law continuum in both the XMM- 
Newton pn and MOS flts and in the PCA + HEXTE -|- 
BAT fits is somewhat low, Fhx ^ 1.5 — 1.6. Using accumu- 
lated RXTE PCA and HEXTE archival data from over 6 
years of monitoring, plus Swift>-BAT 9-month survey data, 
we constrain the strength of the Compton refiection com- 
ponent i? to be < 0.5. We also find the first evidence for 
a high-energy continuum cutoff in this source, at 90±20 
keV. The high-energy continuum is discussed further in 
§9.3. 

The narrow Fe Ka emission line at 6.4 kcV is resolved 
in the pn spectrum. From time-resolved spectral fits to 
the RXTE-PCA data, we find tentative evidence for a sig- 
nificant fraction of the Fe K line flux to track variations 
in the continuum flux F2-10 on time scales of <700 days. 
The 6.4 keV Fe Ka line properties are discussed further in 
§9.4. 

In addition, we find significant evidence for emission 
near 6.0 keV in both the pn and MOS 1-1-2 spectra. It 
is modeled approximately equally well as a narrow emis- 
sion line or as the red wing to a relativistically-broadened 
Fe K emission line, as discussed in §9.5. 

9.1. Overview of the X-ray absorbers 

In the XMM-Newton RGS spectrum, we model two 
zones of absorption, with log^io = l-2lio!o8 ^^'^ log^hi = 
2.90lg'26, with very similar column densities, A^h.Io = 
l.ltg;^ X 1021 cm-2 and iVn.hi = 2.4t2;0 x lO^i cm'^. 
Other AGN have been reported to host multiple warm 
absorber components with different ionization parameters 
but similar hydrogen column densities (see e.g., Blustin et 
al. 2005). The blueshift relative to systemic of the high- 
ionization absorber was significantly constrained, yielding 
an outflow velocity relative to systemic of -(20601^7°) 
s~i. The best-fit blueshift relative to systemic for the low- 
ionization absorber corresponded to an outflow velocity 
relative to systemic of -(420l^gQ) km s~^; however, as this 
is primarily constrained by only two (possibly blended) ab- 
sorption features, models with blucshifts ranging from zero 
to identical to that for the high-ionization absorber were 
not clearly rejected. It is therefore not clear if the two 
zones are physically and kinematically distinct or if there 
exists a single ionized absorber spanning a broad range 
of ionization levels, a possibility to consider if the ionized 
X-ray absorbing gas is spatially extended. For instance, 
Gongalves et al. (2006) model the warm absorbing layers 
in NGC 3783 as a single constant-density gas component 
in pressure equilibrium. 

To estimate the distance r between the central black 
hole and the outflowing gas in NGC 3227, we can use 
^ = Lion/{nr'^), where n is the number density, and Lion 
is the 1-1000 Ryd illuminating continuum luminosity. We 
estimate the maximum possible distance to the material 
by assuming that the thickness Ar must be less than the 
distance rmax- The column density A^h = nAr, yielding 
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the upper limit rmax < iion/(-^HC)- We estimate the un- 
absorbed 1-1000 Ryd flux to be 3 x 10""' erg cm'^ 
s~^, which corresponds to iion ~ 1.5 x iC^ erg s~^ (using 
a distance of 20.3 Mpc, and assuming Ho = 70 km s~^ 
Mpc"i and Aq = 0.73). Ass uming two physically distinct 
warm absorbers, for the high- ionization X-ray absorber, 
this yields rmax = 3.6 pc. For the low-ionization X-ray 
absorber, the constraints are even weaker: rmax = 150 pc. 
These constraints are much poorer compared to that de- 
rived by Gondoin ct al. (2003) for a one-zone absorber, 
''max = 0.45 pc, as those authors found a higher value for 
the column density and a lower value for the ionizing flux. 

However, we can also derive a minimum radial distance 
from the black hole rmin for the winds via the require- 
ment for the outflow velocity v to be greater than the 
escape velocity: For the high-ionization X-ray absorber, 
rmin = {c^/v'^)Rg = 17000i?g = 40 light-days, placing it 
outside both the BLR and the inner radius of the dust as 
reverberation-mapped by Suganuma ct al. (2006). 

Under the assumption that the gas is in equilibrium and 
that the outflow velocity is a constant, the mass outflow 
rate Mout of the X-ray absorbers can be derived via conser- 
vation of mass: Mout = flnr^vnip, where v is the outflow 
velocity, nip is the proton mass, and O is the covering frac- 
tion. We then substitute nr^ = iion/C- Assuming f2 = 0.3 
(arbitrary), we find, for the high-ionization absorber, Mom 
~ 2 X 10^^ gm s~^ ~ 0.03 Mq yr~^. For the low-ionization 
absorber, assuming v = 420 km s~^, Mout ~ 2 x 10^'' gm 
^ 0.3 Mq yr"-'^. We note that the actual outfiow rate 
should be lower if there is an extreme degree of collimation 
along the line of sight. We can compare the outflow rate to 
the inflow accretion rate Mace using Lboi = '?MaccC^, where 
T] is the accretion efficiency parameter, typically 0.1. Woo 
& Urry (2002) estimate L^oi for NGC 3227 to be 7.2 x lO'^^ 
erg s~^, which means an accretion rate relative to Edding- 
ton, Ly,oi/LEdd, of 1.4%. For = 0.1, M^cc = 0.01 M© 
yr""'^. The kinetic power associated with the outflow com- 
ponent, estimated as Mout'y^, is thus in the approximate 
range 10*°~*^ erg s~^. 

The outflow mass rate and kinetic energy are rough es- 
timates only, but it does appear likely that the outflows 
represent at least a large fraction of the AGN's accretion 
rate. If the X-ray ionized outflows are long-lived, then a 
sustained feeding of the black hole would be difficult. 

9.1.1. Helium-like emission line diagnostics 

It is not obvious whether the emission lines detected 
in the RGS spectrum are due to collisional ionization or 
photo-ionization. We do not detect a strong Fe XVII L 
3d-2p -'^Pi line at 826 cV, an indicator of collisional ioniza- 
tion. However, RRC lines, indicators of photo-excitation, 
are not obvious, either. 

We can attempt to use diagnostics associated with the 
helium-like O and N emission triplets. In a collisionally- 
(photo-) ionized plasma, the resonance (forbidden) line 
dominates. However, in NGC 3227, there is not signif- 
icant dominance of one of these lines in either emission 
triplet. The presence of hclium-likc absorption lines may 
be a factor: the emission and absorption lines are not com- 
pletely separated in energy given the modeled blueshift 
of the absorption lines, the RGS resolution and the sig- 
nal/noise ratio of the RGS data. This fact also complicates 



our use of the density indicator R = f /i and the temper- 
ature indicator G = {f + i)/r (e.g., Porquet & Dubau 
2000) where /, i and r are the intensities of the forbidden, 
intercombination, and resonance lines, respectively (z is 
the summed intensity of both intercombination lines). As 
G depends on r, it can indicate whether photo-ionization 
dominates or whether a photo-/collisional-ionization hy- 
brid is applicable (Porquet & Dubau 2000). Porter & Fer- 
land (2007) warn that although G indicates temperature 
in a collisionally-ionized plasma, it should not be used as 
a temperature indicator in a photo-ionized plasma (in ad- 
dition, see Porter & Ferland 2007 for warnings regarding 
usage of R as density indicator in photo- ionized plasmas). 

For the O VII emission triplet in NGC 3227, we mea- 
sure G = 4.7t4 o; G < 2.0 for the N VI triplet. Uncertain- 
ties here are statistical only and do not include systematic 
effects associated with the presence of absorption lines. 
These values thus do not yield any useful constraints on 
whether collisional- or photo-ionization dominates. For 
the O VII emission triplet, R = 1.4^0 g. Assuming the 
plasma is either purely ionized or a hybrid of coUisionally- 
and photo-ionized material, and assuming a temperature 
near 10® K, electron densities of roughly cm~^ 
are implied (Porquet & Dubau 2000). For either collisional 
or photo-ionization, densities above 10^^ cm~^ are ruled 
out by the presence of the strong O VII (/) line. 

9.1.2. Long-term variability in the cold and ionized 
absorbers 

Variations in the column density of neutral absorbing 
gas along the line of sight has been reported to occur on a 
wide range of time scales (days to years) in both Scyfcrt 1 
AGN (e.g., I Zw 1, Gallo ct al. 2007; NGC 4151, Puccetti 
et al. 2007; NGC 3516, Turner et al. 2008, Markowitz et al. 
2008) and Seyfert 2 AGN (Risaliti et al. 2002, 2005). NGC 
3227 is no exception, given the 2000-1 obscuring event by 
a lowly-ionized (log ^ ~ 0) dense cloud of column den- 
sity 2.6 X 10^^ cm~^ and inferred to be located in the BLR 
(Lamer et al. 2003). That cloud is likely physically distinct 
from the local cold gas and the ionized X-ray absorbers, 
as iVn, local, A^H,WAi and AfH,WA2 modeled from the 2006 
XMM-Newton observation are much lower than that for 
the 2000-1 obscuring cloud. Nonetheless, there does seem 
to be evidence for variations in A^h, local on time scales of 
years in NGC 3227: George ct al. (1998b) noted A^H.iocai to 
increase from ^ 3 x 10^° cm"^ in 1993 to ~ 3 x 10^^ cm~^ 
by 1995. Gondoin et al. (2003) reported iVn.iocai = 7x lO^o 
cm~^ from the 2000 XMM-Newton observation; that value 
is the same as the column here in the 2006 EPIC-pn spec- 
trum. Assuming that discrepancies in values of iVn. local 
inferred from different missions are intrinsic to the source 
and not due to cross-instrumental calibration uncertain- 
ties, the range in measured A^h, local values would suggest 
the presence of cold gas along the line of sight at <<pc 
radii from the black hole. 

George et al. (1998b) noted an increase in the ionized 
absorber column density, from ~ 3 x 10^^ cm~^ in 1993 
to ^ 3 X 10^^ cm~^ 1995, ruling out a location outside 
the NLR. The value of A^h.wa obtained by Gondoin et al. 
(2003) using xstar modeling is 2.7±0.7x 10^^ cm"^. Vari- 
ability in Nu,wA between 2000 and 2006 is thus implied to 
be significant only at the 2.3cr confidence level. Further- 



X-ray Spectrum of NGC 3227 



13 



more, Gondoin ct al. (2003) did not include a component 
to model a soft X-ray excess; if the soft excess was in- 
trinsically present during the 2000 observation, then their 
estimate of A^h.wa would be too high. Assuming the vari- 
ations in iVn.wA since 1993 to be real, then at least some 
fraction of the ionized X-ray absorbing gas is likely present 
at <<pc radii. 

9.1.3. Possible connections between the X-ray and UV 
ionized absorbers 

The extended dusty lukewarm absorber (hereafter 
DLWA) studied by Crenshaw et al. (2001) may be a phys- 
ically and/or kinematically distinct component from the 
ionized X-ray absorbers. Crenshaw et al. (2001) report 
the ionization parameter U of the DLWA to be 0.13 (see 
George et al. 1998a for definitions of the optical/UV ion- 
ization parameter U and the X-ray ionization parameter 
C/x). Assuming an optical- to-X-ray spectral index aox 
(estimated from photometric measurements taken from 
NED), and using the conversions from Figure 1 of George 
et al. (1998a), U/U^ ~ 250 and ^/U^ ~ 5000, we can trans- 
late our best-fit ionization parameters for the low- and 
high-ionization X-ray absorbers, log^=1.21 and log^=2.90, 
respectively, into corresponding values of C/ = 1.0 and 40. 
The values of the ionization parameters for the DLWA 
and the low-ionization X-ray absorber are not too dissim- 
ilar, tentatively suggesting a possible physical connection. 
However, it is not possible to definitively link the DLWA 
and the low-ionization X-ray absorber in velocity space, 
as, other than the Fe UTA and O VII edge, there are no 
strong X-ray absorption features >0.35 keV endemic to 
only the low-ionization, and not the high-ionization, X- 
ray absorber at the best-fit ionization parameter values. 

Crenshaw et al. (2001) note that the column density of 
the DLWA is 2 x 10^^ cm~^, very similar to that measured 
for both X-ray absorbers here. The outflow velocities rela- 
tive to systemic of DLWA lines (tens of km s^^) are much 
lower than that of at least the high-ionization X-ray ab- 
sorber (~ 1500 — 2000 km s~^). Furthermore, it is likely 
that the DLWA has a high covering fraction (Crenshaw 
et al. 2001). We cannot rule out the possibility that the 
low-ionization X-ray absorber may be the inner edge of 
the extended DLWA. Moreover, it is conceivable that the 
outflowing, ionized X-ray absorbers may supply the ex- 
tended, 100-pc scale DLWA gas (at least along the line of 
sight), with the velocity slowing and ionization parameter 
decreasing as the distance from the black hole increases. 
If the X-ray ionized absorbing gas is continuously feed- 
ing the UV absorber, we might expect the existence of an 
"intermediate" zone of absorbing material with outflow ve- 
locity between those of the observed X-ray and UV ionized 
absorbers. Confirmation of such a zone would strengthen 
the physical connection between the X-ray and UV ionized 
absorbers. However, clarification of whether or not the X- 
ray absorbing outflow is sustained or intermittent would 
be needed to determine the exact physical and kinematic 
connection between the UV and X-ray absorbers. 

The properties we have derived for the ionized X-ray 
absorbers are consistent with the notion that dust may be 
swept up by an outflowing wind at radii > a few light-days 
(outside the BLR, where dust cannot survive). However, 
the current data cannot distinguish between a dusty and 



a dust-free X-ray ionized absorber; as discussed in §5, the 
upper limit on the neutral Fe L3 dust edge in the RGS 
spectrum corresponds to a hydrogren column density con- 
sistent with that implied by optical line absorption. 

9.2. The variable soft excess in NGC 3227 

The origin of the soft excess emission in NGC 3227 is not 
immediately obvious. The soft excess does not undergo the 

same obsciiration by cold material as that suffered by the 
hard X-ray continuum. However, the large discrepancy 
between Fsx and Fhx argues against both partial cover- 
ing scenarios and against the soft X-rays being nuclear 
power-law emission scattered in an extended region. We 
have modeled the soft excess in NGC 3227 using a black- 
body component to represent direct thermal emission from 
accretion disk, but such a model seems to be unphysical 
given the apparent consistency of blackbody parameters 
across a wide range of Seyfert properties, including black 
hole mass (e.g., Gierlinski & Done 2004). We have also 
modeled the soft excess as inverse Comptonization of seed 
photons, likely thermal optical and UV photons from the 
accretion disk, by thermal electrons (Sunyaev & Titarchuk 
1980). Such a mechanism could conceivably operate in the 
hot, ionized surface of the accretion disk (Hubeny et al. 
2001; Janiuk, Czerny & Madejski 2001). 

In Seyferts, soft excesses have been seen to vary on time 
scales of weeks, e.g., as seen in the narrow line Seyfert 
Is Ark 564 and Ton S180 (Edelson et al. 2002). The 
soft excess in the quasar 3C 273 has also been known 
to vary on time scales of ~ a week (Kim 2001 and ref- 
erences therein). Vaughan et al. (2002) found the soft 
excess of the NLSyl Ton SI 80 to vary rapidly, though in 
concert with the hard X-ray power law component, leading 
to measurents of Fyar roughly independent of energy band 
across the EPIC bandpass. However, most other Seyferts' 
Fvar spectra peak near 1-2 keV (e.g., Arevalo et al. 2008; 
Vaughan & Fabian 2004). This is commonly explained 
with a variable power-law component superimposed over 
a constant or relatively less variable hard component above 
~5 keV (likely the Compton reflection hump) and a con- 
stant or less variable soft component below ~1 keV (the 
soft excess). In NGC 3227, however, _Fvar increases with 
decreasing energy below 1 keV (Figure 14) due to a soft 
excess that is strongly variable in normalization. 

The simultaneous increases in UV and soft X-ray con- 
tinuum flux are not statistically significantly correlated, 
so a direct UV-soft X-ray connection is tenuous at best. 
However, if such a connection were real, and if both the 
UV and soft X-ray originate in the same location on the 
disk, a rapid change in illumination of the accretion disk 
could explain the simultaneous increases. 

We can investigate if the observed UV trend can be at- 
tributed to thermal reprocessing of the increasing soft X- 
ray flux trend by the accretion disk. Given the black hole 
mass, accretion rate, and assuming a standard thin disk 
with an inner radius of 6 i?g, 90% (50%, 10%) of the 260 
nm continuum emission originates from within 60 (25, 10) 
i?g. The light crossing time across the diameter of the disk 
at a 60 Rg radius is roughly 20 ks. In a thermal reprocess- 
ing scenario, the UV flux should track the soft X-ray flux 
and be smeared on a time scale < 20 ks, assuming that 
the soft X-ray continuum emitting region is located near 
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the central disk plane. If the soft X-ray emitting region 
is located well off the plane, the light travel time to the 
UV-emitting part of the disk is increased (e.g., 30 ks travel 
time for a height of 100 i?g.) The light curves displayed 
in Figure 1 arc consistent with this scenario, as any soft 
X-ray to UV lag is less than several tens of ks. 

Another possibility is that the soft X-ray continuum 
flux originates in inverse Comptonization of UV seed pho- 
tons from the disk; e.g., the scenario associated with the 
COMPST model. However, the fact that the observed in- 
crease in the soft X-ray flux (40%) is stronger than the ob- 
served increase in the UV flux (10%) would argue against 
all of the soft X-ray photons being passively-reprocessed 
UV photons, unless only the soft X-ray photons were 
anisotropically beamed along the line of sight, or unless the 
soft X-ray flux is responding to a comparative increase in 
the UV continuum flux that occurred less than a few tens 
of ks before the start of the XMM-Newton observation. 

In the inverse Comptonization scenario, and in the ab- 
sence of thermal reprocessing, the UV variability would 
have to be intrinsic to the disk. We can consider inwardly- 
propagating variations in the local mass accretion rate 
(e.g., Lyubarski 1997, Arcvalo & Uttley 2006), specifically, 
the case where the fluctuations travel on the viscous time 
scale. To produce the 10% observed variation, such a prop- 
agating fluctuation would have to modify 20% of the flux 
contained within a radius of 25 i?g, or 100% of the flux con- 
tained within a radius of 10 i?g. Let us consider the former 
case: a fluctuation that propagates with an inward veloc- 
ity that takes it from a radius of 25 i?g to the innermost 
stable disk radius of 6 i?g over 100 ks, while modulating 
20% of the flux within 25 i?g, could produce the observed 
UV variability. However, the required velocity, ~l/5 of the 
orbital velocity, is very high for a standard disk (Shakura 
& Sunyaev 1973): to be the viscous velocity for a standard 
disk, one would need {H/Kf'a = 0.2 {H, R, and a are the 
height, radius and viscosity parameters, respectively), i.e., 
a very thick disk with a large viscosity parameter, which is 
probably not adequate for the optical-emitting portion of 
the disk. The UV variability is thus likely not associated 
with accretion rate fluctuations traveling on viscous time 
scales. The sound speed near 25 Rg (assuming a = 0.1) 
is roughly two orders of magnitude slower than the or- 
bital speed, and sound waves are thus likely inadequate to 
produce the required flux modification in 100 ks. 

The rapidly variable soft excess is qualitatively similar 
to that observed in 3C 120 with Suzaku by Kataoka et al. 
(2006), who model the variable soft excess in that object as 
the high-energy part of a synchrotron jet component. The 
rapid variability observed in the soft X-ray and UV bands 
in NGC 3227, and the fact that the soft excess imdergoes 
less absorption by cold material compared to the hard X- 
ray continuum, may also be consistent with existence of 
such a jet component, in addition to and independent of 
the accretion disk corona (hard X-ray) emission. Unlike 
3C 120, NGC 3227 is of course a radio-quiet AGN. A -40 
pc scale jet was detected in NGC 3227 by the VLA in 
a 1991 observation (Kukula et al. 1995), though 10-100 
pc scale jets are common in Scyferts (e.g., Gallimore ct 
al. 2006). Broad band spectral energy distribution (SED) 
modeling of an accretion disk -|- jet is beyond the scope 
of this paper; constraining such an SED component would 



be complicated by the presence of absorbing dust in the 
poorly-studied EUV band. However, such a component is 
plausible if its peak in vL^ space were near IQ^'^-^^ Hz 
(explained if the jet were yoimg, e.g., lO'"^"^ yrs), and 
the total observed synchrotron power were no more than 
]^q4i-42 gj.g g-i_ ^]^jg gcgnario, the observed UV contin- 
uum emission would be a mixture of slowly- varying "big 
blue bump" emission from the accretion disk plus rapidly 
variable synchrotron emission (which would contribute less 
than —50% to the total observed UV continuum flux). 

9.3. The high energy emission components 

The photon index of the power-law component used to 
model the hard X-ray continuum emission in both the 
XMM-Newton and RXTE data is rather low (Fhx ~ 
1.5 — 1.6) compared to traditional or "canonical" values 
of - 1.8 - 1.9 for most broad- line Seyfert 1 AGN. How- 
ever, similarly low values of Fhx have been reported pre- 
viously for this object (e.g., George et al. 1998b). The 
presence of a cutoff near 100 keV, as implied by both the 
Swiftr-BAT and the RXTE-REXT'E data, would suggest 
thermal Comptonization, but better data in the 50-200 
keV band and above are needed to confirm the cutoff and 
better constrain its energy. The low-energy rollover to the 
hard X-ray continuum component was modeled above as 
being due to cold absorption, but could also potentially be 
a signature of Comptonization, as it is possible for Comp- 
tonized continuua to have a low-energy rollover e.g., below 
1 keV depending on parameters such as the distribution 
of input photon field (e.g., Titarchuk 1994). 

The 7.1 keV edge is likely due to reflection, as the ab- 
sorbing components modeled here predict an edge depth 
much too small to detect here. However, the strength of 
the Compton reflection component R is measured to be 
low (< 0.5), suggesting that Compton-thick material ex- 
ists in only a small fraction of the sky as seen from the 
continuum source. For example, the putative Compton- 
thick molecular torus invoked in standard Seyfert 1/2 uni- 
fication schemes (Urry & Padovani 1995) could be weak 
or small. Another possibility is that the optically-thick 
accretion disk could be truncated, or the inner accretion 
disk could transition into an optically-thin, radiativcly in- 
efficient flow, e.g., an ADAF or RIAF (Narayan & Yi 
1994, 1995; Blandford & Begelman 1999; Narayan et al. 
2000). Such flows arc frequently invoked when describing 
low accretion rate accreting black hole systems (e.g., low 
luminosity AGN). In addition, Liu et al. (2007) forward a 
model for low accretion rate flows in which optically-thin 
coronal matter condenses into a cool, optically-thick inner 
disk, for LBoi/iEdd = 0.1 — 2% (assuming values for the 
viscosity parameter a — 0.1 — 0.4). Condensation radii of 
a few tens of Rg are plausible. 

Wu & Gu (2008) noted that for both black hole X-ray 
binary systems and AGN accreting above a critical "transi- 
tion" value of iBoi/^Edd, Lhx correlates with Lboi/ L^dd- 
Below this transition, Fhx and iBoi/^Edd anti-correlate. 
Wu & Gu (2008) suggest that the transition between a thin 
disk and an ADAF-typc flow occurs near Leoi/iEdd = 1% 
for black hole X-ray binary systems and near 0.3% for 
AGN. The corresponding value of Fhx at this transi- 
tion point is 1.5, similar to that observed in NGC 3227 
(I/Boi/-^Edd = !%)• This suggests that the accretion flow 
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in NGC 3227 may be consistent with either a thin disk, a 
ADAF-type flow, or a transition to both. 

9.4. The Narrow Fe Ka line at 6.4 keV 

The measured energy of the Fe Ka emission hne is con- 
sistent with an origin in neutral material. The line is re- 
solved in the pn spectrum; the measured width a of 65±14 
eV translates into a FWHM velocity ufwhm of 7000±1500 
km s~^. This velocity width is similar to that for the H/3 
hne, 5500±500 km s'^ (Wandel et al. 1999), suggesting 
an origin for the Fe K line in material commensurate with 
the BLR. 

Assuming Keplerian motion, and assuming the veloc- 
ity dispersion is related to the FWHM velocity as <v^>^ 
I'^^FWHM (Netzer et al. 1990), assuming a black hole mass 
Mbh of 4.22 ± 2.14 x IO'^Mq (Peterson et al. 2004), we 
can use GA/bh = r<v^> to estimate the radius r of the 
line-emitting material. We find r = 7.2;t4'9^ light-days, 
equivalent to 3000t^?^EJ Rg. This value is consistent with 
the radius of < 700 light-days inferred from the very crude 
reverberation mapping discussed in §6.2. It is also consis- 
tent with the BLR (H/3) region radius of 9^g light-days 
(Peterson et al. 2004), as well as with the ~5-20 light- 
day inner radius of the dust as reverberation-mapped by 
Suganuma et al. (2006). 

The Fe line-emitting material may be spatially extended 
in NGC 3227; only the innermost regions have had time to 
respond to continuum variations, which may get smoothed 
out by the outer regions anyway. Alternatively, we note in 
Figure 16 that extrapolation of the -F'2-io~-^FcKa relation 
to zero continuum flux results in a non-zero offset. That is, 
there may exist Fe-line emitting material which does not 
respond to the continuum flux, and only some fraction of 
the total observed Fe line flux responds to the continuum 
variations on < 700 light-days. Either of these two ideas 
may explain why i^var for the 2-10 keV continuum flux is 
larger than that for the Fe line. 

The observed line EW is 91±10 eV. The column density 
of the gas required to produce such an emission line must 
be greater than 10""^^ cm~^, otherwise the optical depth 
would be insufficient to produce such a prominent line. 
However, the lack of an obvious Compton shoulder (the 
emission near 6.0 keV was not identified as such) would 
argue against the bulk of the line photons originating in 
Compton-thick material. 

Let us first assume an origin in optically-thin gas that 
completely surrounds a single, isotropically-cmitting con- 
tinuum source (covering fraction fc — I) and is uniform in 
column density. We can relate EW to Nn using the fol- 
lowing equation, taken from e.g., Markowitz et al. (2007): 

P{E)a^^{E)N^dE 

uo is the fluorescent yield: the value for Fe, 0.34, was taken 
from Kallman et al. (2004). /kq is the fraction of photons 
that go into the Ka line as opposed to the K/3 line; this 
is 0.89 for Fe I. A is the number abundance relative to 
hydrogen; solar abundances, using Lodders (2003), were 
used. P{E) is the spectrum of the illuminating contin- 
uum at energy E] E'linc is the Ka emission line energy. 

^® http:/ /www. pa. uky.edu/~verner/photo. html 



(Tph(-E) is the photo- ionization cross section assuming ab- 
sorption by K-shell electrons only; all cross sections were 
taken from Veigele (197jJ3)- The observed Fe Ka EW can 
thus be produced from material with A^h = 1.4±0.2 x 10^^ 
cm^^. Similarly, the observed Ni Ka line EW of 131|^q eV 
indicates A^h = 3.9t|^ x 10^3 cm'^. 

If the absorbers are instead clumpy and lie out of the 
hne of sight, we can use Eqn. 1 of Wozniak et al. (1998), 
which gives the expected Fe Ka line intensity for a cloud 
with a column > lO^'^ cm^^ lying off the line of sight 
and subtending a fraction VI/Att of the sky as seen from 
a single, isotropically-emitting continuum source. Assum- 
ing solar abundances and given the observed Fe Ka line 
intensity, for values of ^/An — 0.1 (0.3), a column density 
of 1.5 X 10^3 (5 X 10^2) cm-2. 

In cither case, the derived column density is higher than 
that of the hard X-ray absorbing material in the SXPL and 
COMPST model fits to the XMM-Newton EPIC data. It 
is, however, closer to the column density of the obscur- 
ing cloud during the 2000-1 obscuring event (Lamer et al. 
2003). Based on the duration of the obscuring event and 
the implied number density, that cloud was inferred to 
be located in the BLR. The similarities in column density 
and inferred location suggest a connection between the Fe 
line-emitting material and the obscuring cloud of 2000-1. 

9.5. Emission at 6.0 keV 

The Fe K emission profiles in both the pn and MOS 
1-1-2 spectra suggest emission redward of the Fe Ka core. 
However, the present data cannot distinguish between a 
model wherein the emission near 6.0 keV is due to a narrow 
emission line or is instead the red wing of a relativistically 
broadened Fe Ka line with an EW near 80 eV. 

9.5.1. A narrow emission feature? 

We first discuss the emission in the context of modeling 
it as a narrow feature. In the pn spectrum, the line width 
cr was <200 eV; the EW was 21^9^ eV. High-resolution 
spectroscopy has yielded evidence for similar narrow emis- 
sion features in roughly 16 Seyferts to date (see Vaughan 
& Uttley 2008 for a review) . These features are commonly 
interpreted as red- or blue-shifted Fe K features associated 
with " hot-spots," localized areas of the accretion disk illu- 
minated by a localized flare just above the disc (possibly 
due to magnetic reconnection), rather than a central illu- 
minator or an extended corona (e.g., Dovciak et al. 2004, 
Goosmann et al. 2007). 

Vaughan & Uttley (2008) point out that these narrow 
emission features are rarely reported at more than 3.0cr 
confldence. Consideration of selection bias and publica- 
tion bias brings skepticism to the validity of these features 
and raises the possibility that a non-trivial fraction of these 
published lines may be consistent with being due to pho- 
ton noise. In the case of the pn spectrum of NGC 3227, the 
Monte Carlo simulations we performed suggested that the 
line was detected at >99.9% confldence (i.e., <0.1% like- 
lihood that it is due to photon noise). However, we cau- 
tion that the Monte Carlo procedure yields an estimate of 
the detection signiflcancc independently of the existence 
of other high-resolution spectra of NGC 3227 and other 
Seyferts. On the other hand, the marginal detection of the 
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emission feature near 6.1 keV in the MOS 1+2 spectrum, 
simultaneously to the pn data, reduces the likelihood that 
the feature is due to photon noise. In addition, the value 
of |Ax^| associated with modeling the line in the pn spec- 
trum is among the largest reported for narrow emission 
Hues (See Table 1 of Vaughan & Uttley 2008). 

If the feature near 6.0 keV is indeed emission from a " hot 
spot" on the accretion disk, we can estimate its radial dis- 
tance from the black hole. For instance, we can assume 
that the observed redshifting is due solely to Doppler shift- 
ing associated with Keplerian motion of the receding side 
of the disk. We use the 3(t energy range of 5.8 to 6.3 keV, 
corresponding to radial velocities along the line of sight 
of 0.02-0. 09c (assuming an origin in neutral Fe). Assum- 
ing the disk is inchned hy i = 30° (45°) with respect to the 
plane of the sky, using Mbh = 4.2 x lO"^ Mq (Peterson et 
al. 2004), and ignoring gravitational redshifting, this cor- 
responds to Keplerian motion at radii of roughly 12-430 
(25-800) i?sch. 

If the flare is co-rotating with the disk, then a minimum 
radius independent of i is defined by the fact that the line 
energy is consistent with being constant and redward of 
6.4 keV for the 100 ks duration of the observation. Any 
orbital period <100 ks is excluded; otherwise the number 
of photons >6.4 keV would be greater than the number of 
photons <6.4 keV. Radii < 9 i?sch are thus excluded. 

An alternate possibility is that the line-emitting mate- 
rial is in free-fall, with no velocity component transverse 
to our line of sight. Yet another possibility is that the 
observed redshifting is purely due to energy losses of pho- 
tons escaping from the vicinity of the black hole; assuming 
for simplicity emission from neutral Fe originating in gas 
lying along the line of sight to the black hole, with all ve- 
locity transverse to line of sight, radii of ~ 6 — 30i?sch are 
plausible. 

Another speculative possibility is that the line-emitting 
material could be associated with the base of the jet, form- 
ing above the accretion disk, as opposed to being a part 
of the accretion disk. Close to the rotation axis in AGN, 
magnetic fields are likely twisted by the differential rota- 
tion of the disk, serving to launch and coUimate jets (e.g., 
Blandford & Payne 1982). For instance, Marscher et al. 
(2008) presented observational evidence for the jet in the 
radio-loud AGN BL Lac to follow a spiral flow as it is 
accelerated through a zone containing a helical magnetic 
field. The 6.0 keV line in NGC 3227 may be emission 
from Fe in a blob of material, with a column density A^h 



Evidence for relativistically broadened Fe line " diskline" 
profiles in a significant fraction of Seyferts has been ac- 
cumulating since the days of ASCA (Tanaka et al. 1995, 
Fabian et al. 2000). In recent years, XMM-Newton and 
Chandra-BETGS observations have demonstrated that a 
narrow core at 6.4 keV is ubiquitous; accurately modeling 
that component, as well as absorption due to ionized gas 
in the line of sight, is critical for accurate measurement of 
the fiux and profile of any broad Fe line present. 

In NGC 3227, our best-fit "DL" model incorporates a 
narrow Fe Ka core, and suggests a relatively weak broad 
line, with an EW of 8II3Q eV. We constrained the in- 
ner radius to be < 22i?g and the inclination to be <25° 
(though these values were derived with the radial emissiv- 
ity parameter fixed at -3). 

Assuming solar abundances, the abundances of Lodders 
(2003), and given the measured strength of the Compton 
reflection component R — 0.40 ± 0.07, the predicted Fe 
fine EW (George & Fabian 19910 is 54±9 eV, consistent 
with the modeled EW of the observed broad Fe line. The 
low value of the broad line EW may also suggest that the 
optically-thick disk is not spatially extended, again consis- 
tent with the notion that the nucleus of NGC 3227 could 
harbor an ADAF or RIAF, and the upper limit on the in- 
ner radius of 22 Rg is potentially consistent with the model 
of Liu et al. (2007) containing a small, optically-thin disk. 

10. CONCLUSIONS 

We have observed the nucleus of the Seyfert 1.5 AGN 
NGC 3227 with XMM-Newton in December 2006 for al- 
most 100 ks. We have used EPIC and RGS spectra to 
study the Fe K bandpass, ionized absorbers, and 0.2-10 
keV continuum in detail. We also present X-ray contin- 
uum light curves and a UV continuum light curve obtained 
with the Optical Monitor. We have combined these data 
with archival RXTE-PCA and HEXTE monitoring, plus 
Swift-JiAT monitoring data, to constrain the level Comp- 
ton reflection, study the high energy continuum up to 200 
keV, and track the flux behavior of the Fe K line on time 
scales of weeks to years. Our main results are summarized 
as follows: 

The EPIC pn spectrum shows the prominent narrow Fe 
Ka emission line to be consistent with an origin in neu- 
tral Fe (Gaussian energy centroid 6.403 ± 0.009 keV). Its 
intensity and EW are 3.5 ± 0.4 x 10~^ ph cm~^ s~^ and 



of 10' 



-21-22 



cm , which has been caught up in the launch- 



91±10 eV, respectively, consistent with an origin in ma- 
terial with a column density io~22.5-23 ^j^-s^ rpj^^ ^^^^^ 
is resolved in the pn spectrum, with a FWHM velocity 
7000±1500 km s"""^. Assuming purely Keplerian motion, 
we estimate the radius of the line-emitting material to be 
7.2j:4^g^ light-days. The FWHM velocity and estimated 
radius are consistent with the BLR (as mapped by H/3 
emission) as well as with the ~5-20 light-day inner radius 
of the dust as reverberation-mapped by Suganuma et al. 
(2006). 

Time-resolved spectral fitting to the 1999-2005 RXTE- 
PCA monitoring data reveal tentative evidence for a signif- 
icant fraction of the Fe K line photons to track variations 
seen in the continuum, with a light travel time delay which 
is not tightly constrained, but is < 700 days. We thus rule 

George & Fabian 1991 predict R = EW/150 eV, using the abundances of Morrison & McCammon (1983), wliicli assumed 3.3 X 10~^ Fe 
atoms per H atom. Lodders (2003) sets the Fe abundance to 2.95 X 10~^ Fe atoms per H atom. 



ing and formation of the weak jet. If the material is close 
to the axis, above the disk, and moving with a velocity 
with a large azimuthal component, then an observed red- 
shift is possible unless the rotation axis lies very close to 
the line of sight. However, the weakness of Seyferts' jets 
may possibly be attributed to weaker acceleration and/or 
collimation mechanisms compared to those in blazars, re- 
ducing the likelihood that such azimuthal velocities can 
exist in Seyferts. 



9.5.2. A relativistically broadened diskline component? 
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out the possibility that the bulk of the (variable) Fe line 
photons originate at distances of 1-2 pc and more, e.g., in 
a multi-pc-scale molecular torus. However, more intensive 
sampling on time scales of days to weeks is required to 
beter constrain the lower limit on the continuum line lag. 

Emission near 6.0 keV is detected in both the pn and 
MOS 1+2 spectra. The fact that it is detected in both 
instruments and the relatively large associated value of 
Ax^ both argue against this feature being an artifact due 
to photon noise. It is inconsistent with being a Comp- 
ton shoulder to the 6.4 keV Fe Ka line. We modeled the 
emission using a narrow Gaussian component, with en- 
ergy ccntroid (in the pn spectrum) 6.04tQ Q4 kcV, width 
a < 200 eV, intensity 9^1 x lO'^ ph cm s'^ and EW of 
2ltl^ eV. A possible origin is a " hot-spot" in the accretion 
disk, originating at radii of tens to hundreds of Rg. The 
emission feature is modeled equally well as the red wing 
to a relativistically broadened "diskline" component, with 
inner radius <22 Rg, inclination <25°, and EW = SltH 
eV. 

Broadband (0.2-10 kcV) EPIC spectral modeling re- 
veals a strong soft excess dominating below 1 keV. In the 
pn spectrum, the soft excess is fit well by a steep power 
law, with Fcjx ^ 3. A blackbody component and a low- 
temperature Cornptonization component each fit the data 
well, though the physical plausibility these latter two com- 
ponents in Seyferts is not generally accepted. The normal- 
ization of the soft excess increases by about 20% during 
the first ^^20 ks of the XMM-Newton observation, and by 
~40% over ~ 50 ks. The soft X-ray band is more strongly 
variable than the hard X-ray band on time scales of tens 
of ks: unlike the Fyar spectra of many other Seyferts, the 
Fvar spectrum of NGC 3227 below 1-2 keV continues to 
increase as photon energy decreases. Such relatively rapid 
variability in the soft excess, independent of rapid vari- 
ability in the hard X-ray band is very rare for Seyferts. 

The OM shows the UV continuum flux (near 260 nm) 
to increase by 10% during the XMM-Newton observation. 
Such <1 day variability is relatively strong compared to 
similar-duration UV continuum light curves obtained with 
the OM (Smith & Vaughan 2007). 

We present the highest-quality gratings spectrum ob- 
tained for NGC 3227 to date. In the RGS spectrum, we 
model the ionized absorption using two zones. The higher- 
ionization zone has log^hi = 2.90lo!26' ouflow ve- 

locity relative to systemic of -(2060l^7o) ^™ '^^^^ 
indicates a minimum radial distance from the black hole 
of 40 light-days, placing it outside the BLR radius. The 
lower-ionization zone has log^io = 1-21^0.08 (corresponding 
to a factor of roughly 3-10 higher than the dusty lukewarm 



absorber modeled by Kraemer et al. 2000 and Crenshaw 
ct al. 2001). Its main signature is a Fe M shell UTA near 
740-780 eV. The best estimate for its outflow velocity rela- 
tive to systemic is -(4201^9°) s~-^. We find no evidence 
for non-solar abundances. Both absorbing zones have col- 
umn densities A^^h.wa near 1 — 2 x 10^^ cm~^. iVn.wA and 
the column densities of local, cold absorption measured 
from the EPIC spectra are too low to be directly associ- 
ated with the obscuring BLR cloud in 2000-1 studied by 
Lamer ct al. (2003). However A'h.wa is similar to that of 
the 100-pc scale DLWA, consistent with the notion that 
the outflowing X-ray absorbing material may supply the 
UV-absorbing DLWA material. 

In the RGS spectrum, we detect five narrow emission 
lines detected at the systemic rcdshift; we identify these 
lines as being due to N VI (r), N VII , and the (f), (i) 
and (r) lines of O VIII. No emission due to Fe XVII L 
line (3d-2p) at 826 eV (a signature of collisionally-ionized 
plasma) is detected (upper limit of 1 eV) . No strong RRC 
features due to H- or He-like ions (signatures of photo- 
ionized plasma) are detected (upper limits range from 5 
to 50 eV). 

We presented the total spectrum derived from RXTE 
PCA and HEXTE archival monitoring data taken in 1996 
and from 1999-2005, and combined them with a 4-channel 
Swift-BAT spectrum from the 9-month survey data (ob- 
tained in 2005). The HEXTE and BAT spectral data re- 
ported here represent the first detailed spectrum of NGC 
3227 above ~ 20 keV and up to almost 200 keV. The hard 
X-ray continuum photon index in both the RXTE+^PlT 
spectra and XMM-Newton EPIC spectra is rather fiat, 
Fhx = 1-5 — 1.6. The strength of the Compton refiec- 
tion hump is rather low, R < 0.5, consistent with the EW 
of the diskline component modeled in the pn spectrum 
(assuming solar abundances). We also find evidence for a 
high-energy continuum cutoff at 90±20 keV. 

The low values of Fhx and R, combined with the low 
value of iBoi/-^Edd ~ 1%, are consistent with the notion 
that NGC 3227 may harbor an optically-thin, radiatively- 
inefficient flow such as an ADAF, in addition to or instead 
of a standard geometrically-thin, radiatively-efflcient disk. 

A.M. thanks M. Elvis, A. Marscher and D. Evans for 
helpful suggestions. This work is based on an observation 
obtained with XMM-Newton, an ESA science mission, and 
has made use of HEASARC online services, supported by 
NASA/GSFC, the NASA/IPAC Extragalactic Database, 
operated by JPL/California Institute of Technology un- 
der contract with NASA, and the NIST Atomic Spectra 
Database. 
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Table 1 

Fe K bandpass components in the time-averaged epic pn spectrum 





P a T a m <3 1 <3 v 

X Oil CLlllCLid 


Value 






983.2/990 


Fe Ka emission line 


Energy (keV) 


6.403 ± 0.009 




Width a (eV) 


65 ± 14 




Intensity (ph cm" ^ s" ^) 


3.5 ±0.4 X 10-5 




EW (eV) 


91 ± 10 


Ni Ka emission line 


Energy (keV) 


~ 41 +0.08 




Intensity (ph cm"^ s~^) 


4+3 X 10-8 




EW (cV) 


,0+13 


6.04 keV emission line 


Energy (kcV) 


fi 04+0-20 




Width cr (eV) 


< 200 




Intensity (ph cm"^ s""*^) 


9.1 ± 2.5 X 10-8 




EW (cV) 


22 ± 6 


Fe K Edge 


Energy (keV) 


7.11 (fixed) 




Optical Depth r 


0.05 ± 0.03 



Note. — Parameters for the best-fitting "GA" model (wherein the emission near 6.0 keV is 
modeled with a narrow Gaussian component), fit to >4 keV data. The model includes a Fe K/3 
line, with energy fixed at 7.056 keV (rest frame), width cr tied to that of the Ka line, and intensity 
equal to 0.13 times that of the Kq line. The width a for the Ni Ka line was also tied to that for 
the Fe Ka line. Equivalent widths EW were determined relative to a locally-fit continuum. 



Table 2 

Model fits to the 0.2-10 keV time-averaged EPIC pn spectrum 



Component 


Parameter 


SXPL Model 


BB Model 


COMPST Model 




xVdof 


1866.9/1724 


1870.4/1724 


1863.4/1723 


Local cold absorption 


A^H.local (cm"2-) 


S.7t°oi X 1020 


5.1+1-2 ^ -^020 


7.7+°-^ X 1020 


Hard X-ray power-law 


Thx 


1.57 ± 0.02 


1 61+0-01 


1 65+''-02 




Norm. (1 keV) 


6.7+°'3 X 10-3 


7.4 ±0.1 X 10-3 


8.1+05 X 10-3 


Hard X-ray Absorption 


JVh.hx (cm-2) 


2.9lo';8 X 1021 


< 6.1 X 1020 


4.2+0-3 ^ ^q21 


Low-ioniz. X-ray Warm Abs. 


Wh.Io (cm-2) 


1.0+Q 3 X 1021 


7.912-7 X 1020 


9.3+2-1 ^ -^q20 




log^io (erg cm s-l) 




1.38 ± 0.19 


1 43+0-20 
l-^'^-0.43 


High-ioniz. X-ray Warm Abs. 


A^H.hi (cm-2) 


1.8+Q g X 102i 


1.8+1-3 X 1021 


1.5+J-3 X 1021 




log^hi (erg cm s-l) 


2 93+0-15 

^•"■^-0.09 


2.9ll°;J? 


2 89+0-1* 


Soft X-ray power-law 


Tsx 


o oc-l-0.27 
•J'-JO-CIO 








Norm. (1 keV) 


4.0l°:2 X 10-3 


83li 

1.7210-37 X 10-4 




Blackbody 


Temperature ksT (eV) 
Norm. 






Comptonized component 


Temperature fesT (keV) 
Optical Depth r 






0.35lH^ 
2412 




Norm. 




58+0-02 


4.4ll0;22 X 10-3 


O VII emission line 


Energy (keV) 


0.58 ± 0.01 


0.58 ±0.01 




Intensity (ph cm-^s 


2.410-5 X 10-" 


1.2l0;8 X 10-4 


2.1 ±0.4 X 10-4 



Note. — SXPL, BB and COMPST refer to the models wherein the soft excess is modeled with a steep power-law component, 
a blackbody component, and a Comptonized component (CompST), respectively. The units of normalization of the power-law 
components arc ph keV^i cm^^ g-i at 1 keV. The blackbody normalization is L-^ijD^g, where L39 is the source luminosity 
in units of lO'*'' erg s^i and Diq is the distance to the source in units of 10 kpc. See the XSPEC user manual for units of the 
CompST component normalization. The reader is reminded that the BB and COMPST models are not generally accepted as 
physically plausible descriptions of soft excesses in Seyferts. The uncertainties listed on the intensity of the O VII emission line 
arc statistical only and do not include likely systematic uncertainties due to blending with narrow absorption lines. The width 
a of the O VII emission line was kept fixed at 0.5 eV. 
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Table 3 

Emission lines fit to the RGS spectrum 



Line 


Energy 


Intensity 


F-test 


Identification 


(cV) 


(ph cm^'^ s^^) 


Ax^ Probability 



N VI (r) 


431 ± 1 


9.9l:fl X 10-5 


-12.0 


0.042 


N VII 


499+2 


8.2t|-^ X 10-5 


-10.4 


0.064 


O VII (f) 




i.33t;;;3,3 X 10-4 


-25.2 


1.4 X 10-3 


O VII (-i.) 


567 ± 1 


9.51?,:? X 10-5 


18.5 


7.8 X 10-3 


O Ml (r) 


o72 ] 


i.oi^-.y'' X 11) " 


7.9 


0.12 



Note. — Best-fit line emission parameters using Gaussian components. All line widths a were 
kept fixed at 0.5 eV. The uncertainties listed on the intensities (Col. [3]) arc statistical only and 
do not include likely systematic uncertainties due to blending with narrow absorption lines. The 
F-test probability values P (Col. [5]) denote the changes that the null hypothesis (line not included 
in model) is valid; the line detection significance is 1 - P. 



Table 4 



Best-fit model parameters for the RGS spectrum 



Component 


Parameter 


Value 




xVdof 


785.7/417 


Local cold absorption 


A^H.local (cm-2) 


10.5l^-g X 10^0 


Soft X-ray power-law 


Tsx 


3.00'±0.25 




Norm. (1 keV) 


6.2+°-3 X 10-3 


Low-ioniz. X-ray Warm Abs. 


■^Vh.Io (cm-2) 


l.ltti X 1021 




logfto (erg cm s-l) 


1 oi-fO.lS 
-^•^-^-0.08 




2lo (absolute) 


+0 00246+'' """4 




Zlo (rel. to systemic) 


-0 00140+° "i'*4 


High-ioniz. X-ray Warm Abs. 


A^H.hi (cm-2) 


2.4+2-0 X 1021 




log^hi (erg cm s-l) 


2 90+0'2i 




Zhi (absolute) 


-0.00302l0;™;!57 




Zhi (rel. to systemic) 


-0.0068810:^1^ 



Note. — Results are shown for the best-fit soft X-ray power-law model, including five emission 
lines due to He-like O and H- and He-like N; see Table 3 for emission line parameters. The units of 
normalization of the power-law components are ph keV-l cm— 2 s— l at 1 keV. See text for details 
regarding uncertainty range in z. The model also included a hard X-ray power- law component 
absorbed by a column ATh^hx ~ 8 x lO^i cm-2 and emerging only above ~ 1.5 keV. 
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Fig. 1. — The top three panels show the EPIC-pii count rate hght curves for the 0.2—10, 0.2—1, and 3—10 keV bandpasses, binned to 600 i 
The bottom panel shows the OM UVWl light curve, binned to 1400 s. 
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Fig. 2. — Data/model residuals to spectral fits to the Fe K bandpass of the EPIC-pn spectrum. Data are rebinned by a factor of 12. 
Vertical dashed lines denote energies of 6.04, 6.40, 7.06, and 7.37 keV. Panel a) shows residuals to a simple power-law model. In panel b), the 
Fe Ka and Fe K/3 emission lines have been modeled. The 7.11 keV edge and Ni Ka line have been modeled in panels c) and d), respectively. 
In panel e), a, Gaussian with an energy centroid of 6.24 keV has been added. However, as shown in panel f), the best-fit "GA" model, with a 
Gaussian with an energy centroid at 6.04 keV, does a superior job in modeling the remaining ~6 keV residuals. The small dip near 6.8 keV is 
not significant and is likely an artifact of fitting. In panel g), the best-fit "DL" model, the 6.04 keV line has been removed and replaced with 
relativistically broadened Fe Ka and Fe K/3 diskline profiles; the narrow Gaussian at 6.04 keV (panel f)) seems to model the residuals better. 
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Fig. 3. — Contour plots showing the results of applying a "sliding Gaussian" (with width a fixed at 10 cV) to a model to the EPlC-pn data 
consisting of a power-law plus a narrow Gaussian at 6.40 keV to model Fe Ka emission. 1-, 2-, 3- and 4-a confidence levels for two interesting 
parameters are denoted by solid, dashed, dot-dashed, and dotted lines, respectively. Residuals near 6.0 keV are clear; additional residuals at 
7.0 and 7.4 keV (above the Fe K edge at 7.11 keV) are investigated further in §3. 
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Fig. 4. — Contour plot of intensity versus line centroid energy for the Gaussian used to model the narrow 6.0 keV emission feature in the 
pn spectrum. The width cr was left as a free parameter. 1-, 2-, 3- and 4-cr confidence levels for two interesting parameters are denoted by 
solid, dashed, dot-dashed, and dotted lines, respectively. 
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Fig. 5. — Data/model residuals to various model fits to the MOS 1+2 spectrum. Data are rebinned by a factor of 8. Panel a) shows 
residuals to a simple power-law model. In panel b), the Fe Kq and Fe K/3 emission lines have been modeled. In panels c) and d), respectively, 
the 6.11 keV emission line and the Fe K edge have been added. 
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Fig. 7. — Contour plot of intensity versus line centroid energy for the Gaussian used to model the narrow 6.1 keV emission feature in the 
MOS 1+2 spectrum. The width a was left as a free parameter. 1-, 2-, 3- and 4-(T confidence levels for two interesting parameters are denoted 
by solid, dashed, dot-dashed, and dotted lines, respectively. 
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Fig. 8. — Data/model residuals to spectral fits to the 0.2-10 keV EPIC-pn spectrum. All models include narrow Gaussian components 
to model the Fe Ka, Fe K/3, Ni Ka and 6.04 keV emission lines, an Fe K edge at 7.11 keV, and a column of neutral gas to model Galactic 
absorption. Panel a) shows the results when a simple hard X-ray power-law component (HXPL) is used. In Panel b), a neutral absorbing 
column has been added to the HXPL, and a steep, soft X-ray power-law component (SXPL) has been added. In panels c) and d), the low- 
and hi-ionization warm absorbers, respectively, have been been included. In panel e), a narrow Gaussian to model O VII emission has been 
included; this is the best-fitting "SXPL" model. In panel f), the soft power-law has been replaced with a blackbody component; this is the 
best-fitting "BB" model. Panel g) shows the results from modeling the soft excess as blurred refiection from an ionized disk. (The small dip 
near 6.8 keV is not significant and is likely an artifact of fitting.) 
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Fig. 9. — Unfolded spectra for the best-fitting SXPL, BB and COMPST models. Panel d) shows the spectra for three selected time-resolved 
segments, further illustrating the increase in soft X-ray flux between segments 1, 2 and 6. 
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Fig. 10. — The RGS spectrum for NGC 3227. The solid line shows the best-fitting unfolded model, described in §5. Data have been 
rebinned by a factor of 2, and are plotted in the rest frame. Black and red data points denote RGS 1 and 2, respectively. The vertical lines, 
which indicate identified lines and edges, denote rest-frame (systemic) energies. The residuals near 687 and 850 eV are instrumental artifacts. 
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Fig. 11.— The upper panel shows the RXTE-PCA (black), RXTE-HEXTE cluster A (red) and B (purple), and Swift-BAT data (blue). 
The middle and lower panels show x residuals to the best-fit models with and without a high-energy cutoff in the power-law component, 
respectively. In each case, the model contains a Compton reflection component and Fe Ka emission. 
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Fig. 12. — Contour plots showing the strength of the Compton reflection component, K (measured via a pexrav component; top panel) 
and cutofi' energy (bottom) versus hard X-ray photon index from spectral fits to RXTE PCA and HEXTE and Swift-BAT data. 1-, 2-, 3- 
and 4-(T confidence levels for two interesting parameters are denoted by solid, dashed, dot-dashed, and dotted lines, respectively. 
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Fig. 13. — The results of fitting the SXPL model to the ten time-resolved spectral segments. From top to bottom: Normalization of the 
soft X-ray power-law component, Fgx, normalization of the hard X-ray power-law component, and Phx- 
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Fig. 14 . — i^var Spectra derived from EPIC-pn data for the entire duration (black), the first 20 ks (red) and the final 80 ks (blue). 
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Fic;. 15. — The results of time-resolved spectral fits to the medium-term (top) and long-term fbottom) RXTE data, showing 2—10 keV 
continuum flux in units of 10~^^ erg cm~^ s~^, Fe K line intensity /peKa ™ units of 10~^ ph cm~ s~^, and Fhx- 
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Fig. 16.- Zero-lag correlation diagrams in whieh Fe line intensity /pcKai in units of 10^^ ph cm~^ s""', is plotted against 2—10 keV 



continuum flux in units of 10 



erg cm 



^. Dashed lines indicate the best-fittine; linear model. 
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Fig. 17. — ICF (red dashed line) and DCF (solid black line) functions for the 0.2-1 keV soft X-ray light curve versus the OM UV continuum 
light curve (positive lag indicates soft X-ray leading the UV.) 



